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1.  SUMMARY 


The  research  completed  under  the  contract  F44620-75-C-0064 
"Research  in  Seismology"  during  the  period  1  July  1978  - 
31  December  1978  falls  within  the  two  broad  topics  of  (1)  Seismic 
source  mechanisms  and  (2)  Seismic  wave  propagation.  The 
specific  research  tasks  within  each  category  have  broad 
applicability  to  the  problem  of  discrimination  of  earthquakes 
from  underground  nuclear  explosions,  with  particular  emphasis 
on  the  Asian  continent. 

A  major  work  on  both  source  mechanisms  and  surface  wave 
propagation  was  Patton's  (1978)  thesis  research  on  the 
simultaneous  determination  of  source  moment  tensor  and  propa¬ 
gation  characteristics  (phase  velocity,  Q)  from  a  collection 
of  sources  in  a  small  region  observed  by  a  set  of  common 
stations.  Taylor  and  Toksoz  (1978)  applied  the  three- 
dimensional  inversion  technique  of  Aki  and  others  to  P 
wave  data  from  the  New  England  seismic  network,  demonstrating 
that  the  distinct  velocity  structure  of  major  geologic  units 
in  an  ancient  continental  collision  zone  extends  to  depths  of 
200  km.  Johnston  (1978)  and  Johnston  and  Toksoz  (1979a, b) 
present  laboratory  data  on  attenuation  in  rocks  at  ultrasonic 
frequencies  to  investigate  several  possible  loss  mechanisms 
for  the  crust,  including  crack  friction  and  fluid  saturation. 

Details  of  these  studies  are  given  in  the  preprints  and 
abstracts  in  the  following  sections.  A  cumulative  list  of 
publications  completed  under  the  contract  is  also  included. 
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2.  SEISMIC  SOURCE  MECHANISMS 


SOURCE  AND  PROPAGATION  EFFECTS  OF  RAYLEIGH  WAVES 
FROM  CENTRAL  ASIAN  EARTHQUAKES 
BY  HOWARD  PATTON 


Submitted  to 

the  Department  of  Earth  and  Planetary  Sciences 
on  April  28,  1978 

in  partial  fulfillment  of  the  requirements  for  the 
degree  of  Doctor  of  Philosophy 


ABSTRACT 

A  reference  point  equalization  method  has  been  developed 
which  enables 'the  separation  of  source  and  propagation  effects  of 
surface  waves.  The  method  works  on  seismic  events  located  in  a 
small  source  region,  which  allows  us  to  assume  that  all  events 
share  the  same  path  effects  to  a  given  receiver.  Two  important 
steps  in  the  method  are  initialization  and  iteration.  Initiali¬ 
zation  obtains  the  first  "reference  events"  in  order  to  compute 
initial  estimates  of  phase  velocity  and  attenuation  coefficient. 
Iteration  simultaneously  refines  the  propagation  parameters  and 
determines  the  source  parameters  of  new  earthquakes  in  the  vicin¬ 
ity  of  the  reference  point.  This  method  was  applied  to  earth¬ 
quakes  in  the  Pamir  mountains,  Central  Asia  (reference  point: 
39.58N,  73.55E). 

In  the  initialization  step,  the  method  of  Weidner  and  Aki 
(1973)  was  applied  to  obtain  focal  depths  and  revise  fault  plane 
parameters  of  the  first  two  earthquakes.  The  residuals  obtained 
from  fitting  the  observed  amplitude  ratios  and  phase  differences 
indicate  that  the  crust  and  upper  mantle  in  Central  Asia  is  more 
laterally  heterogeneous  than  near  the  ocean  rifts,  the  site  of 
Weidner  and  Aki's  experiment.  We  computed  heterogenity  quotients 
of  .46  vs  .80  x  lO-^Napier^/km  and  .93  vs  2.4  x  10-4radian2/km 
for  ocean  versus  continent  as  a  measure  of  the  increasing  scatter 
in  amplitude  and  phase  of  20- 60sec  Rayleigh  waves  due  to  lateral 
heterogeneities. 

To  determine  source  parameters  in  the  iteration,  we  applied 
the  linear  moment  tensor  inversion  on  Rayleigh  wave  complex 
source  spectra.  The  presence  of  random  additive  errors  in  the 
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complex  spectra  does  not  pose  difficulties  for  recovering  reli¬ 
able  source  parameters  using  the  linear  inversion  method.  Howev¬ 
er,  amplitude  magnification  errors  in  the  complex  spectra  will 
lead  to  over-estimation  of  the  moment  tensor  elements  and  phase 
incoherency  will  lead  to  under-estimation.  In  applying  this 
method  to  our  dataset,  it  was  necessary  to  modify  the  straight 
least  squares  inversion  method  because  of  its  sensitivity  to  even 
a  few  bad  data  points. 

The  residuals  obtained  from  the  repeated  application  of  the 
moment  tensor  inversion  over  trial  focal  depth  showed  two  minima: 
one  minima  at  depths  less  than  20km  and  the  other  at  depths 
greater  than  70kra.  The'values  of  the  residuals  at  these  minima 
were  close  enough  to  cast  doubt  on  the  determination  of  focal 
depth.  One  way  to  r.esolve  this  ambiguity  is  to  compare  the 
geometry  of  the  moment  tensor  obtained  for  shallow  and  deep  focus 
inversions  with  observed  P-wave  polarities. 

The  focal  depths  of  eight  out  of  the  nine  events  in  our  da¬ 
taset  were  found  to  be  shallow,  between  5-15km.  Their  principal 
compressive  stress  axes  are  aligned  north-south  and  nearly  hor¬ 
izontal,  consistent  with  the  interpretation  of  plate  tectonics  in 
Central  Asia.  The  moment  tensor  inversion  generally  gives 
three-couple  force  systems  having  significant  non-zero  intermedi¬ 
ate  component.  However,  in  light  of  errors  in  our  data  and  the 
similarities  of  the  double  couple  models  to  the  three-couple 
models,  it  can  not  be  established  convincingly  that  these  results 
are  caused  by  departures  of  source  from  the  double  couple  model. 

We  interpret  the  propagation  parameters  in  terms  of  lateral 
variation  of  phase  velocity  and  Q  on  the  Eurasian  continent.  A 
phase  velocity  regionalization  is  proposed  involving  five  con¬ 
tinental  provinces:  Indian  Shield,  Northern  Platforms,  Coastal 
Plains,  Tectonic  and  Plateau.  Phase  velocities  on  the  Indian 
Shield  are  20%  higher  than  velocities  on  the  Plateau  province  at 
26sec  period  and  5%  higher  at  90seo  period.  Stable  provinces  in 
Eurasia  are  found  to  have  significantly  higher  phase  velocity 
than  tectonic  provinces  out  to  150sec  period. 

Interpretations  of  the  phase  velocities  on  the  Indian  Shield 
show  a  lithosphere  thickness  of  about  120km,  considerably  thicker 
than  the  lithosphere,  under  the  Northern  Platforms  (*v75kra) .  The 
lowest  shear  velocity  in  the  upper  mantle  is  found  under  the  Tec¬ 
tonic  province  with  a  value  about  4.3km/sec  over  depths  between 
83-240km.  The  crustal  thickness  of  the  Plateau  province  is  as 
great  as  70km  provided  that  shear  velocity  in  the  lower  half  of 
the  crust  is  about  3.8km/sec.  The  upper  mantle  structure  under 
this  province  is  very  similar  to  James’  (1971)  for  the  Andes 
mountains  region. 

Surface  wave  amplitudes  on  the  Eurasian  continent  are 
strongly  affected  by  horizontal  refraction  as  well  as  intrinsic  Q 
of  the  medium.  Average  Rayleigh  wave  Q  has  very  different  char¬ 
acter  on  paths  over  northern  platforms  and  tectonic  provinces 
east  and  west  of  the  reference  point.  Under  platforms  Q  is  found 
to  Increase  with  depth  from  a  Q  between  200-300  in  the  crust  to 


300-500  in  the  lid  and  asthenosphere.  West  of  the  reference 
point  Q  is  low  (^60)  in  the  bottom  20km  of  the  crust  and  extend¬ 
ing  n.35km  into  the  mantle.  East  of  the  reference  point  Q  appears 
to  be  %60  in  a  layer  35km  thick  at  the  base  of  the  crust.  Deeper 
in  the  mantle  Q  increases  suggesting  that  the  low  Q  zone  is  con¬ 
centrated  at  shallow  depths  in  the  upper  mantle  under  the  tecton¬ 
ic  provinces  of  Asia. 
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Pip.3-l.lJ  Azimuthal  equidistance  projection  of  Eurasia  centered  on 
the  reference  point  in  the  Pamir  Mountains.  All  station 
used  in  this  study  are  shown  on  this  map. 
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3.  SEISMIC  WAVE  PROPAGATION 


Three-Dimensional  Crust  and  Upper  Mantle  Structure 
of  the  Northeastern  United  States 


Steven  R.  Taylor  and  M.  Nafi  Toksflz 

Department  of  Earth  and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


November,  1978 


Three-Dimensional  Crust  and  Upper  Mantle  Structure  of  the 

Northeastern  United  States 

Abstract 

Teleseismic  P  wave  arrival  times  recorded  by  the  North- 
Eastern  Seismic  Network  are  used  to  invert  for  lateral  crust 
and  upper  mantle  structure  to  depths  of  500  km.  Three-dimen¬ 
sional  inversion  of  the  travel  time  data  between  two  ancient 
erogenic  provinces  suggests  that  structures  down  to  possibly 
200  km  and  greater  can  be  correlated  with  surficial  geologic 
and  tectonic  features.  This  has  the  important  implication 
that  major  orogenic  belts  have  effects  that  reach  well  into 
the  lithosphere  which  are  stable  for  extended  periods  of  time 
perhaps  as  long  as  1  billion  years.  The  crust  beneath  the 
Paleozoic  Appalachian  Province  is  characterized  by  slightly 
greater  thicknesses  and  lower  average  velocities  than  that 
of  the  Precambrian  Grenville  Province;  The  higher  average 

velocities  associated  with  the  Grenville  province  extend  to 

♦ 

depths  of  200  km  and  appear  to  be  maximum  beneath  the  Adi¬ 
rondack  dome.  A  relatively  low  velocity  anomaly  extending 
to  depths  in  excess  of  200  km  and  dipping  to  the  northwest 
shows  a  spatial  correlation  with  the  Bronson  Hill  -  Boundary 
Mountains  Anticlinorium  in  central  New  Hampshire  and  Maine. 
These  structures  occupy  the  sites  of  a  complex  series  of 
island  arc  sequences  last  active  in  Early  Devonian  time  prior 
to  the  Acadian  orogeny.  This  low  velocity  region  may  repre¬ 
sent  subducted  oceanic  lithosphere  which  has  undergone  post- 
orogenic  radioactive  heating. 
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Introduction 

This  paper  reports  on  the  preliminary  results  of  a  study 
of  teleseismic  P  wave  arrival  times  recorded  by  the  Northeastern 
Seismic  Network.  The  inversion  of  teleseismic  P  wave  travel 
time  residuals  is  a  commonly  used  method  of  extracting  velocity 
information  in  the  crust  and  upper  mantle  beneath  am  array. 
Notable  studies  include  those  in  California  (Bolt  and  Nuttli, 
1966;  Nuttli  and  Bolt,  1969;  Steeples  and  Tyer,  1976)  at 
LASA  (Sheppard,  1966;  Iyer  and  Healy,  1972;  Greenfield  and 
Sheppard,  1964) ,  in  Arizona  (Johnson,  1967) ,  in  Japan  (Zandt, 
1975;  Hirahara,  1977) ,  at  NORSAR  (Aki,  1977) ,  in  Yellowstone 
National  Park  (Iyer,  1975;  Zandt,  1978),  in  Hawaii  (Ellsworth 
and  Koyanagi,  1977),  in  New  Madrid  (Mitchell  et  al. ,  1977), 
and  beneath  the  Tabela  Array  (Menke,  1977) . 

Many  of  the  earlier  studies  consisted  of  determination 
of  residual  patterns  and  development  of  qualitative  models 
to  explain  the  observed  trends.  Recent  models  have  utilized 
the  techniques  developed  by  Aki  et  al.  (1977) 7  for  mapping 
the  three-dimensional  structure  beneath  an  array. 

To  date,  very  little  seismological  work  pertaining  to 
crust  and  upper  mantle  structure  has  been  done  in  the  north¬ 
eastern  United  States.  Leet  (1941) ,  Lin eh am  (1962) ,  and  Katz 
(1955)  published  refraction  results  using  a  small  permanent 
array  combined  with  some  portable  stations.  More  recently 
some  refraction  work  has  been  reported  by  Chiburis  and  Graham 
(1978)  and  Aggarwal  in  Schnerk  et  al.  (1976) .  Other  relevant 
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Geologic  Setting  and  Tectonic  Evolution 

The  northern  Appalachians  have  undergone  a  long  and 
extremely  complex  geological  development ,  Dating  and  strati¬ 
graphic  correlation  between  regions  has  been  complicated  by 
numerous  orogenic  episodes .  The  northern  Appalachians  can 
be  divided  into  three  major  tectonic  units  (Bird  and  Dewey, 
1970;  Naylor,  1975) :  a  western  belt  and  an  eastern  belt 
possible  representing  the  margins  of  two  once  convergent 
continental  masses  surround  a  central  orogenic  belt  composed 
mainly  of  eugeoclinal  lithologies  (Figure  1) .  The  western 
unit  is  mainly  underlain  by  rocks  of  the  Precambrian  Grenville 
Province  which  are  exposed  in  the  Adirondacks  and  outlying 
massifs  such  as  the  Green  Mountains,  Berkshires,  and  the 
Hudson  Highlands .  Unconformably  overlying  the  Grenville 
basement  is  an  Eocambrian  to  Cambrian  platform  sequence  which 
grades  upward  into  a  Lower  Ordovician  clastic  sequence.  Found 
above  are  the  Taconic  klippen  which  consist  primarily  of 
deep-water  shales,  sandstones,  and  graywackes.  Paleontological 
evidence  suggests  that  they  were  deposited  contemporaneously 
with  the  shelf  sediments  (Zen,  1972). 

The  central  orogenic  belt  consists  of  a  number  of  broad 
structural  warps.  The  Connecticut  Valley  Synclinorium  (CVS) 
exists  to  the  east  of  the  previously  discussed  Precambrian 
massifs  and  can  be  traced  from  Connecticut  through  Quebec  to 
the  Gulf  of  St.  Lawrence.  The  CVS  contains  a  thick  highly 
metamorphosed  eugeoclinal  sequence  divided  into  two  members 


separated  by  a  major  Middle  Ordovician  unconformity.  A 
linear  serpentinite  belt  follows  the  western  flanks  of  the 
CVS. 

East  of  the  CVS  lies  the  Bronson  Hill  Anticlinorium  (BHA) 
which  consists  of  a  chain  of  elliptical  gneissic  domes 
(Oliveran  Plutonic  series  in  New  Hampshire) .  The  structure 
can  be  traced  from  Connecticut  through  northern  New  Hampshire 
and  is  probably  continuous  with  the  Boundary  Mountains  Anti¬ 
clinorium  in  Maine.  Mantling  the  domes  is  a  series  of  mafic 
metavolcanics  associated  with  felsic  metavolcanics  and  meta¬ 
sediments  (Ammonoosic  Volcanics  in  New  Hampshire)  of  Middle 
Ordovician  age  or  older  (440±30  m.y.;  Naylor,  1975).  Uncon- 
formably  overlying  the  Ammonoosics  is  a  series  of  Silurian 
to  Lower  Devonian  highly  metamorphosed  elastics  with  some 
carbonates  and  volcanics. 

Eastward  of  the  BHA  lies  the  Merrimack  Synclinorium  (MS) 
which  is  a  major  northeast-trending  tectonic  feature  extending 
from  eastern  Connecticut  through  Maine  and  into  New  Brunswick. 
It  is  the  site  of  thick  accumulations  of  Ordovician  to  Lower 
Devonian  metasediments  typically  metamorphosed  to  sillimanite 
grade.  These  metasediments  (Devonian  Littleton  Formation)  can 
be  correlated  with  Devonian  strata  at  the  top  of  the  BHA.  The 
MS  also  contains  large  volumes  of  intrusives  belonging  to  the 
Middle  Devonian  New  Hampshire  Plutonic  Series  and  the  Mesozoic 
White  Mountain  magma  series. 

On  the  eastern  flank  of  the  MS  is  a  major  northeast¬ 
trending  thrust  belt  (Clinton-Newbury,  Bloody  Bluff,  and  Lake 
Char  Faults)  extending  from  southern  Connecticut  through 
eastern  Massachusetts  (Skehan,  1969).  Magnetic  anomalies 


associated  with  the  formations  in  the  thrust  belt  suggest 
that  the  faults  continue  offshore  in  an  east-northeast 
direction  into  the  Gulf  of  Maine  (Weston  Geophysical,  1976; 
Alvord  et  al. ,  1976)  and  possibly  into  New  Brunswick  (Nelson, 
1976) .  In  eastern  Massachusetts  the  northwest  dipping 
thrusts  greatly  offset  metamorphic  isograds  and  no  strati¬ 
graphic  units  can  be  traced  across  them  (Pat  Barosh,  personal 
commun. ) .  Deformation  and  metamorphism  in  the  MS  probably 
delineates  a  zone  of  mavimtim  intensity  of  the  Middle  Devonian 
Acadian  Orogeny  (Rodgers,  1970). 

The  Eastern  Basement  is  exposed  to  the  east  of  the  above 
described  thrust  belt.  These  units  are  probably  correlative 
with  rocks  of  the  Avalon  Zone  in  Newfoundland  and  southeastern 
New  Brunswick  (Bird  and  Dewey,  1970;  Nelson,  1976).  The 
region  in  eastern  Massachusetts  is  characterized  by  plutonic, 
metasedimentary  and  mstavolcanic  rocks  metamorphosed  mainly 
to  chlorite  grade  ranging  in  age  from  late  Precambrian  to 
Carboniferous.  Onfossiliferous  strata  and  scattered  age 
dating  has  made  geologic  interpretations  enigmatic  (Naylor, 

1975;  Zarrow  et  al. ,  1978).  However,  according  to  Naylor 
(1975) ,  no  rocks  have  been  assigned  an  age  greater  than  650  m.y. 
which  is  significantly  younger  than  the  Grenville  age  rocks 
(M.,100  m.y.)  in  the  western  belt. 

A  possible  scenario  for  the  evolution  of  the  northern 
Appalachians  will  be  outlined  below.  It  must  be  emphasized  that 
this  is  an  over-simplified  first  order  model,  based  on  often 
poorly  defined,  conflicting  geology.  However,  many  of  the 


large  scale  features  are  correlative  with  observed  velocity 
anomalies  described  in  a  later  section. 

The  oldest  rocks  in  the  study  area  are  exposed  in  the 
Adirondacks  and  the  Precambrian  outliers  of  the  Green  Moun¬ 
tains,  Berkshires,  and  Hudson  Highlands.  On  the  basis  of 
petrological  and  structural  arguments  Dewey  and  Burke 
(1973)  suggested  that  the  Adirondacks  represent  deep  erosional 
levels  of  part  of  an  ancient  Tibetan-type  plateau  formed 
behind  a  continental,  collision  zone  during  the  Grenville 
Orogeny  (^1,100  m.y.).  Doming  in  the  southern  Adirondacks 
may  be  due  to  intersection  of  regional  fold  structures  affecting 
Precambrian  rocks  (McLelland,  1977) .  Since  the  late  Precambrian 
the  Adirondacks  have  apparently  remained  a  regional  high,  as 
suggested  by  the  onlap  of  Paleozoic  sediments  (Rodgers,  1970) . 

It  is  currently  thought  that  continental  rifting  stage 
initiated  approximately  820  m.y.  ago  leading  to  the  formation 
of  the  Iapetus  Ocean  (proto- Atlantic)  (Rankin,  1976) .  The 
late  Precambrian  and  Cambrian  geology  of  the  western  belt 
is  characterized  by  the  establishment  of  an  Atlantic-type 
stable  continental  margin. 

Late  Precambrian  lithologies  on  the  eastern  belt  also 
indicate  a  rifting  stage  with  the  development  of  an  active 
continental  margin  (Kennedy,  1976).  Geochemical,  paleo- 
magnetic,  and  paleontological  evidence  suggests  that  the 
western  and  eastern  belts  were  located  on  opposite  sides 
of  the  Iapetus  Oceam  (Strong  et  al . ,  1974;  McKerrow  and 
Cocks,  1976;  Kent  and  Opdyke,  1978). 

Early  or  Mid-Ordoviciam  through  Permiam  times  are 


characterized  by  the  episodic  closing  of  the  Iapetus  Occam. 

The  BEA  was  a  site  of  major  volcanic  activity  in  this  time 
period  as  evidenced  by  the  presence  of  thick  volcamic  sequences. 
Major  deformation  occurring  in  Middle  and  Late  Ordoviciam 
time  marked  the  climax  of  the  Taconic  Orogeny  which  affected 
rocks  within  amd  west  of  the  BHA.  The  Taconic  klippe  were 
thrust  at  this  time  and  later  deep  seated  thrusts  involved 
Grenville  basement  resulting  in  the  emplacement  of  the  Pre- 
cambrian  massifs.  The  linear  belt  of  ultramafics  found  to 
the  east  of  the  massifs  were  emplaced  at  this  time,  and 
probably  represent  highly  altered,  Cambrian  to  Lower  Ordovician 
(Chidester,  1968)  obducted  oceanic  crust  of  the  Iapetus  Ocean 
or  of  a  marginal  sea  behind  the  volcanic  arcs  of  the  central 
mobile  belt.  The  sequence  of  tectonic  and  metamorphic  events 
are  indicative  of  effects  due  to  long-term  heating,  increased 
ductility  and  basement  reactivation  due  to  induced  convection 
above  subducted  lithosphere.  This  effect  is  observed  in 
continent-continent  convergence  zones  (Toksdz  and  Bird,  1975) 
as  well  as  on  active  continental  margins  involving  only 
subductive  oceanic  lithosphere  and  arc-continent  collisions 
(Burchfiel  and  Davis,  1975). 

Evidence  based  on  styles  of  deformation  and  me tamorphl sm 
indicate  that  the  Taconic  Orogeny  in  New  England  probably  was 
an  episode  of  arc-continent  collision.  Numerous  island  arc 
segments  and  inner-arc  basins  now  located  within  the  central 
mobile  belt  were  probably  involved  in  the  deformation.  The 
lack  of  any  major  unconformity  and  Ordovician  deformation  in 
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the  MS  (Moench  and  Zartman,  1976)  indicates  that  oceanic 
terrain  continued  to  separate  continental  crust  of  the 
eastern  belt  from  North  America.  However,  the  amount  and 
type  of  volcanism  in  the  central  mobile  belt  following  the 
Taconic  Orogeny  indicate  that  closure  of  the  Iapetus  was 
nearly  complete.  It  may  have  been  the  close  approach  of 
the  eastern  belt  and  possibly  the  attempted  subduction  of 
an  oceanic  ridge  which  caused  the  island  arcs  to  converge 
onto  North  America. 

The  Silurian  to  Early  Devonian  time  was  a  period  of 
relative  quiescence  and  erosion  of  the  Taconic  highlands 
(Boucot,  1968).  In  the  Early  Devonian  era  there  is  evidence 
of  increased  tectonic  activity  as  indicated  by  the  deposition 
of  vast  thicknesses  of  turbidite  sequences  across  the  major 
synclinoria  (Littleton  Formation  of  New  Hampshire)  and 
renewed  volcanism  along  arcs  in  the  central  mobile  belt. 

This  episode  of  major  deformation  is  called  the  Acadian 
Orogeny  and  appears  to  have  spanned  a  period  of  approximately 
-30  m.y.  (Naylor,  1971)  which  is  similar  to  the  duration 
of  the  main  collisional  stages  of  the  Himalayas  (Toksdz  and 
Bird,  1977). 

The  Acadian  orogeny  is  typified  by  the  following  over¬ 
lapping  sequence  of  events:  1)  high  grade  metamorphism  in 
southeastern  New  England,  2)  intrusion  of  granitic  plutons, 

3)  large  scale  westward  recumbent  folding  of  many  of  the 
plutons,  and  4)  brittle  deformation  and  development  of  large 
thrust  belts.  These  styles  of  deformation  show  many  parallels 
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with  those  documented  in  other  continental  convergence  zones 
(Dewey,  1977).  The  period  of  metamorphism,  intrusion,  and 
ductile  deformation  are  deep  crustal  processes  characteristic 
of  a  major  thermal  event  associated  with  initial  heating  of 
the  crust  from  below  due  to  asthenospheric  convection  prior 
to  collision  and  post-collisional  radioactive  heat  generation 
(Toksdz  and  Bird,  1977) .  The  binary  granites  of  the  New 
Hampshire  Plutonic  Series  (dated  around  360  m.y.;  Naylor, 

1975)  sure  indicative  of  extensive  crustal  heating  and  partial 
fusion  of  crustal  sediments  and  are  similar  to  those  found 
throughout  much  of  the  Himalayas  (Bird,  1976) . 

Three  major  events  following  the  Acadian  Orogeny  are 
recorded  in  New  England.  Alleghenian  deformation  in  Pennsyl¬ 
vanian  to  Permian  times  strongly  affected  rocks  in  southeastern 
New  England.  The  Allegheny  Orogeny  probably  represents  the 
final  closure  of  the  Xapetus  Ocean  by  the  collision  of 
Africa  with  North  America.  It  has  been  suggested  (Dewey  and 
Kidd,  1974)  that  the  Alleghenian  deformational  belt  truncates 
the  Appalachian-Caledonian  belt  in  southern  Connecticut  and 
Rhode  Island  and  is  correlative  with  the  Hercynian  belt 
in  central  Europe. 

“  In  late  Triassic  time  large  rift-valley  graben  systems 
were  developed  along  the  trend  of  the  Appalachians .  These 
grabens  were  filled  with  thick  accumulations  of  non-marine 
red  beds.  The  rifting  was  accompanied  by  the  fissure  eruption 
of  basaltic  lavas  and  the  intrusion  of  numerous  diabase 
dikes,  sills  and  stocks.  This  period  of  time  represents  the 
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initial  opening  of  the  modern  Atlantic  Ocean.  The  most  recent 
major  tectonic  activity  in  New  England  was  the  emplacement 
of  the  White  Mountain  Magma  Series,  intruded  over  a  100  m.y. 
interval  from  Jurassic  to  Cretaceous  time  (Chapman,  1976). 

The  northeastern  U.S.  is  presently  an  intraplate  region 
adjacent  to  a  stable  continental  margin,  characterized  by 
relatively  low-level,  diffuse  seismicity.  It  is  not  clear 
at  this  time  what  structures  control  seismicity  patterns. 

The  region  within  and  possibly  east  of  the  Grenville  Province 
is  characterized  by  a  northeast  trending  maximum  compressive 
stress  and  some  earthquakes  appear  to  be  controlled  by  re¬ 
activated  paleozoic  faults  (Sbar  and  Sykes,  1977) .  However, 
it  is  not  known  whether  this  stress  pattern  extends  into 
New  England. 


•  v  v  v  v  .  v 


222 


’>  \»  ->  v  ■>  ’.l  •>w.i  v  v  v  v  •>  ■j'">  v  v  w  ?.v  v  .■  -.-.-I 


22 


Travel  Time  Residuals 


Teleseisaic  P-wave  arrival  times  were  collected  from  50 
stations  which  are  part  of  the  Northeastern  Seismic  Network 
(Fig.  2) .  -The  stations  were  installed  and  are  operated  by 
research  groups  at  Boston  College  (Weston  Observatory) , 

Lament- Doherty  Geological  Observatory,  and  Massachusetts 
Institute  of  Technology.  Station  distribution  is  fairly 
uniform,  although  there  are  areas  with  large  station  spacing. 
Not  all  of  the  stations  were  in  operation  for  the  duration 
of  this  study.  Because  the  network  was  not  in  full  operation 
until  February,  197*,  events  used  in  this  study  covered  a 
time  range  between  February,  1976,  and  December,  1977. 

Most  stations  are  equipped  with  vertical  component  short 
period  1  or  2  Hz  geophones.  Signals  are  telemetered  to 
central  recording  sites  (located  at  each  of  the  above  insti¬ 
tutions)  and  recorded  on  develocorder  film. 

The  relative  arrival  times  of  teleseismic  P-waves  were 
read  from  enlarged  copies  of  16mm  develocorder  film.  In 
general,  strong  scattering  in  the  frequency  bandwidth  of 
the  incident  P  waves  is  not  observed  as  indicated  by  the 
coherency  of  the  first  few  cycles  across  the  network.  The 
signals  from  each  of  the  stations  were  visually  correlated 
using  an  overlay  of  representative  waveforms.  Relative  arrival 
measurements  were  taken  from  a  prominent  peak  or  trough  early 
in  the  signal,  and  readings  were  taken  to  the  nearest  0.1  sec. 

A  misidentification  of  the  peak  or  trough  would  cause  a 
residual  error  of  0.5  to  1.0  sec.  which  is  readily  apparent 
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once  the  residuals  are  calculated. 

1822  arrival  times  from  68  events  were  read  from  stations 
of  the  Northeast  network.  The  epicentral  distances  ranged 
between  25  0  and  95  °  and  azimuths  were  mainly  from  the  north¬ 
east  (from  Eurasia) ,  south  and  southwest  (Central  and  South 
America),  and  the  northwest  (Aleutians,  Kuriles),  (Figure  3). 

Elevation  corrections  were  applied  to  the  data  and  travel 
times  were  redtrced  to  a  datum  elevation  at  sea-level.  Absolute 
P-wave  residuals  were  calculated  using  Herrin  travel-time 
tables  and  are  defined  to  be 


R?j  -  T!js  -  Ti?  «» 

TJ 

where  R^  is  the  absolute  residual  with  respect  to  Herrin 
tables  for  station  i,  event  i;  T?.  is  the  observed  travel 
time  using  event  locations  and  origin  times  from  PDE  bulletins; 

n 

is  the  theoretical  travel  time  through  a  Herrin  earth. 

The  residuals  were  further  reduced  by  calculating  relative 
residuals  with  respect  to  a  mean  residual  computed  for  each 


event ; 


R.  .  *  R?.  -  J  I  R?. 
X]  13  N  13 


where  N  is  the  number  of  stations  reporting  P  arrivals  for 
event  j. 

Using  the  above  relationships,  positive  residuals  repre¬ 
sent  late  arrivals  where  the  incident  rays  have  been  slowed 
in  the  crust  or  uppermantle  beneath  the  network.  It  should  be 
noted  that  no  significant  differences  were  observed  in  the 
values  of  the  relative  residuals  when  Jeffreys-Bullen  travel 


time  tables  were  used  instead  of  Herrin  Tables.  An  error 
in  the  reported  origin  time  or  a  source  correction  term 
will  affect  all  stations  similarly  and  will  drop  out  when 
relative  residuals  are  computed.  The  effects  of  the  structure 
at  the  source/  heterogeneities  encountered  along  the  major 
portion  of  the  travel  path,  and  structural  variations  beneath 
the  array  are  major  causes  of  the  absolute  residuals. 
Intuitively,  the  cone  of  rays  from  a  teleseismic  source  to 
different  stations  in  the  network  will  sample  nearly  identical 
travel  paths  until  they  diverge  in  the  upper  mantle  and  crust 
beneath  the  array.  Thus,  it  is  reasonable  to  assume  that 
the  computed  average  residuals  contain  information  regarding 
this  common  travel  path  up  to  the  point  where  the  rays  diverge 
significantly.  By  subtracting  off  the  average  residual, 
effects  due  to  mislocation  errors  or  heterogeneities  encountered 
along  the  common  path  are  greatly  reduced. 

The  above  approximations  which  are  reasonable  when  network 
aperture  is  small  become  progressively  less  valid  as  the 
network  aperture  is  increased.  If  the  apparent  network 
aperture  is  large  enough,  the  curvature  (second  derivative) 
of  the  travel  time  curve  becomes  significant.  The  problems 
involved  in  the  analysis  of  relative  P-wave  residuals  have 
been  discussed  by  Engdahl  et  al.  (1977).  Using  data  from 
Alaskan-Aleutian  stations,  Engdahl  concluded  that  the  primary 
source  of  observed  variations  in  relative  residuals  was  from 
local  structure  beneath  the  network.  For  station  separations 
comparable  to  the  maximum  diameter  of  the  Northeastern  U.S. 
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Network  (600  km) ,  scatter  in  relative  residuals  due  to 
epicenter  mislocations ,  source  structure,  earth  models,  etc. 
were  at  most0.2  sec,  which  is  little  more  than  our  reading 
accuracy.  Analysis  of  relative  residuals  recorded  in  the 
northeastern  U.S.  supports  -these  conclusions.  We  found  no 
systematic  differences  in  relative  residuals  from  various 
depths  in  the  same  epicentral  regions  where  it  is  presumed 
ra^s  would  encounter  different  source  structures. 

Two  additional  sources  of  error  remain  to  be  discussed. 
Because  of  down-time  for  some  stations,  and  of  other  stations 
that  became  operational  during  the  period  of  this  study, 
slightly  different  station  subsets  are  available  for  each 
earthquake,  and  the  mean  residual  subtracted  off  from  a 
given  source  region  may  vary  for  different  events.  However, 
this  problem  is  greatly  reduced  by  the  large  number  of 
stations  (typically  25-35)  used  to  calculated  the  mean  residual 
for  each  earthquake .  The  other  problem  results  from  the  fact 
that  each  subnetwork  records  on  a  different  time  base.  Each 
sub-network  operates  a.  temperature  compensated  crystal 
oscillator  clock  and  timing  corrections  are  generally  less  than 
20  ms/day.  To  assure  accuracy  tests  were  made  by  comparing 
signals  from  a  common  station  recorded  on  more  than  one  sub¬ 
network.  The  analysis  of  travel  time  residuals  were  carried 
out  using  .relative  residuals.  The  distribution  of  average 
relative  residuals, R,  for  each  station  is  shown  in  Figure  4. 

The  relative  residuals  as  a  function  of  azimuth  and  incidence 
angle  show  significant  variations  across  the  network  (Figure  5) . 


These  trends  suggest  the  presence  of  large  scale  lateral 
heterogeneities  in  the  crust  and  upper  mantle,  and  considerable 
effort  is  required  to  separate  the  factors  causing  the  ob¬ 
served  patterns. 

Average  relative  residuals  for  each  station  were  calculated 
with  ten  or  more  readings  from  well-distributed  azimuths. 

Large  positive  residuals  (late  arrivals)  occur  throughout 
central  and  northern  New  Hampshire,  southern  Maine  and 
eastern  Vermont.  As  was  observed  by  Fletcher  et  al. (1978) 
these  contrast  sharply  with  large  negative  residuals  in 
western  Vermont,  northern  New  York,  southeastern  New  York, 
and  southeastern  Connecticut.  Shallow,  localized  structural 
differences  would  result  in  rapidly  varying  trends  between 
adjacent  stations.  However,  the  slowly  varying  distribution 
of  residuals  suggests  that  the  observed  variations  are 
probably  due  to  deep,  regional  structures  such  as  differences 
in  crustal  thickness  and  upper  mantle  velocity. 

As  a  first  step  we  compare  the  residuals  with  surface 
geologic  and  other  geophysical  data.  Figure  7  illustrates 
the  relationship  between  average  relative  station  residuals, 

K,  and  the  regional  Bouguer  gravity  taken  from  each  station 
site  (Kane  et  al.  (1972),  Simmons,  (1964),  and 

Diment  et  al.  (1972) ) .  Because  the  long  wavelength  (slOO  km) 
regional  gravity  anomalies  probably  reflect  variations  in 
Moho  topography,  (Simmons,  1964),  the  inverse  correlation 
between  gravity  and  average  residuals  suggests  that  for  many 
stations,  the  average  residuals  depend  partly  on  crustal 
thickness  variations.  It  is  also  interesting  to  note  the 
dependence  of  average  residuals  on  the  age  of  the  basement 
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beneath  each  station.  The  two  sets  of  symbols  in  the  figure 
correspond  to  the  inferred  age  of  the  basement  for  each 
station  site.  Stations  in  New  York  and  New  Jersey  are 
located  above  Precambrian  (>1.Q  b.y.)  Grenville  basement, 
while  most  New  England  stations  are  atop  younger  (<1.Q  b.y.) 
Paleozoic  basement  ‘of  the  Appalachian  province  (the  central 
mobile  belt,  and  eastern  belt  described  earlier) .  Clearly, 
the  stations  within  the  older  Grenville  Province  see  early 
arrivals  relative  to  those  in  the  younger  Appalachian  Province. 
However,  because  of  the  relationship  between  Bouguer  gravity 
anomalies  and  station  elevation  it  is  difficult  to  attribute 
this  separation  totally  to  crustal  thickness  variations 
(Figure  8).  In  isostatically  compensated  regions,  crustal 
thickness  should  be  correlated  with  surface  elevations 
(Woolard,  1954;  Tseng,  1975).  Figure  8  illustrates  that 
there  exists  a  correlation  between  gravity  and  elevation 
although  there  is  no  obvious  separation  between  the  two  regions. 
To  explain  the  behavior  of  both  Figures  7  and  8  it  is  necessary 
to  invoke  velocity  changes  in  addition  to  any  possible 
variations  in  crustal  thickness.  This  idea  is  supported  by 
refraction  profiles  constructed  in  the  northeastern  U.S. 

Although  little  detailed  refraction  work  has  been  done  in 
this  region,  existing  studies  indicate  slightly  higher  average 
crustal  velocities  in  the  Grenville  Province,  as  opposed  to 
the  Appalachian  Province.  Table  I  shows  MIT  and  Weston 
Geophysical  Co.  (Ron  Street,  personal  comm.)  refraction 
models  collected  from  measured  quarry  blasts  for  eastern 
Massachusetts,  southern  and  central  New  Hampshire ,  and 


southern  Maine.  Also  shown  is  a  model  published  by  Aggarwal 
in  Schnerk  et  al.  (1976)  for  northern  New  York  and  the 
Adirondacks.  Comparison  of  the  refraction  models  indicates 
that  the  crust  beneath  this  part  of  the  Appalachian  Province 
is  on  the  average  slightly  thicker  and  slower,  and  it  is 
possible  to  account  for  a  large  part  of  the  almost  one  second 
difference  in  relative  residuals  between  the  two  regions. 

The  dependence  of  average  relative  residuals  on  crustal 
properties  is  also  supported  by  the  results  of  a  time-term 
analysis  which  indicated  a  positive  correlation  between  Pn 
residuals  and  teleseismic  P-wave  residuals  (Taylor  and  ToksOz , 
1978) .  However,  these  observations  cannot  fully  explain 
Figures  7  and  8  and  as  indicated  by  three-dimensional  models 
presented  in  the  next  section,  structural  differences  between 
the  two  regions  extend  into  the  upper  mantle. 

In  addition  to  distinct  patterns  in  the  distribution  of 

•  • 

average  residuals  across  network,  most  stations  show  a  strong 
dependence  in  residuals  with  azimuth  and  incidence  angle. 

Polar  plots  of  relative  residuals  as  a  function  of  azimuth 
and  incidence  angle  for  many  northeastern  stations  are  shown 
in  Figure  5.  Moving  from  south  to  north  we  see  a  number  of 
interesting  residual  trends.  In  Connecticut,  for  arrivals 
from  the  northeast,  the  residuals  become  increasingly  less 
positive  (arrive  earlier)  as  the  incidence  angle  decreases 
until  in  southwest  Connecticut  and  northern  New  Jersey  residuals 
are  largely  negative.  Rapidly  varying  structure  is  observed 
to  the  south  which  will  be  reflected  in  the  three-dimensional 
models  presented  in  the  next  section.  In  Massachusetts , 
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arrivals  are  late  to  the  northwest  and  small  to  the  south  and 
northeast  while  New  Hampshire  and  southern  Maine  observe  late 
arrivals  from  all  azimuths.  Negative  residuals  are  observed 
in  the  Grenville  Province  across  most  of  New  York,  partic¬ 
ularly  in  the  northern  regions,  where  large  negative  residuals 
are  observed  from  the  northwest. 

The  residual  patterns  indicate  that  between  the  Grenville 
ana  Appalachian,  orogenic  belts  structural  differences  exist 
in  the  crust,  and  into  the  upper  mantle.  It  appears  that 
a  large  region  of  relatively  low  velocities  occurs  in  the 
upper  mantle  beneath  central  new  England  and  that  velocities 
become  increasingly  higher  beneath  the  older  Grenville  Province 
to  the  west.  From  Figure  4  it  is  evident  that  most  stations 
in  close  proximity  to  one  another  show  a  similar  distribution 
in  residual  values  implying  that  lateral  variations  are  net 
a  clear  surface  effect.  Since  the  incidence  angles  are  small 
for  the  teleseismic  P-waves  (15  # -  30°  )  ,  lateral  variations 
in  crustal  velocities  cannot  account  for  the  observed  effects. 
For  example,  teleseismic  rays  of  equal  slowness  and  incident 
from  opposite  azimuths,  enter  the  crust  with  a  maximum  sepa¬ 
ration  of  50  km.  Azimuthal  variations  of  one  second  and 
greater  would  require  unreasonably  large  lateral  velocity 
differences. 

The  possibility  of  a  dipping  Moho  as  a  cause  of  azi¬ 
muthal  variations  in  residuals  was  investigated.  To  the 
first  order  in  the  dip,  a  ,  the  maximum  differential  travel 
time  for  teleseismic  ray  approaching  from  opposite  azimuths 
perpendicular  to  the  strike  of  a  dipping  layer  is  given  by 
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where  h  is  the  depth  to  the  interface  beneath  a  given 
station;  a  is  the  dip  of  the  interface;  9  *  sin"1  (pVm) 
represents  the  incident  angle  at  the  interface;  and  n  * 

w 

^  .  Substitution  into  equation  (3)  implies  that  unreasonably 

c 

large  regional  dips  of  the  Moho  (20-3(f  )  are  required  to 
satisfy  the  observations.  Thus,  it  is  reasonable  to  assume 
that  laterally  varying  structures  in  the  upper  mantle  be- 
neatHthe  network  are  required  to  explain  the  pronounced 
residual  versus  azimuth  and  incident  angle  variations. 


Three-Dimensional  Inversion 


Qualitative  description  of  residual  patterns  as  in  the 
previous  section  is  a  useful  but  inadequate  means  of  esti¬ 
mating  the  structural  disturbances  an  incident  ray  encounters 
beneath  the  network.  In  this  section  we  present  results 
from  a  three-dimensional  inversion  of  the  teleseismic  travel 
time  data  following  the  techniques  illustrated  by  Aki,  et  al 
(1977).  The  initial  model  used  in  the  inversion  consists 
of  a  plane-layered  medium  where  each  layer  is  sub-divided 
into  a  series  of  right  rectangular  blocks.  Using  relative 
teleseismic  P  wave  residuals,  individual  slowness  (^r)  per¬ 
turbations  are  computed  for  each  block  relative  to  a  layer 
average.  Aki  et  al.  (1977)  has  shown  that  for  the  ith 
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station  and  jth  event,  relative  residuals  (with  respect  to 
a  mean)  can  be  represented  by 
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,  reduced  travel  times 


-  travel  time  through  kth  block 
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D, ..  -  -  r j-,  percent  perturbation  in  slowness 
13  vk  relative  to  the  layer  average  (1/V^) 

N  -  number  of  observations  for  event  j 

a  -  number  of  layers 


In  matrix  form 
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where  r  -  n- length  vector  containing  relative  residuals 

m  -  k- length  vector  containing  unknown  velocity 
perturbations 

A  -  n  x  k  symmetric  matrix  with  reduced  travel  times 

Solution  of  the  system  (5)  by. classical  least  squares 

fails  because  there  exist  two  sources  of  non-uniqueness. 

The  first  source  can  be  caused  by  insufficient  data  and 

results  in  vertical  "coupling"  of  inadequately  sampled  blocks 

in  separate  layers  which  share  a  common  ray  path.  The  other 

source  of  non-uniqueness  is  inherent  to  the  formulation  of 

the  problem  and  results  from  the  removal  of  the  average 

residual  from  the  data.  This  implies  that  within  a  given 

layer,  the  velocity  perturbations  are  determinant  only  to 

within  a  constant  and  no  information  regarding  absolute 

velocities  is  attained.  Thus,  the  system  (5)  is  solved  by 

using  a  damped  least  squares  technique  where  a  positive 

T 

constant  is  added  to  the  diagonal  elements  of  A  A  which  reduces 
the  effect  of  zero-eigenvalues  on  the  generalized  inverse 
solution..  The  solution  to  (5)  is  given  by 

fi  -  <ATA  +  02!)"1  ATr  (6) 

where  0*  is  the  damping  parameter  and  Z  is  the  identity  matrix. 
The  resolution  and  covariance  matrices  are,  respectively 

R  *  <ATa  +  0 *1) "1  ATA  (7) 

C  -  S*  (ATA  +  02!)"1  R  (8) 

where  S2  is  the  estimate  of  the  data  variance,  o^2. 


Because 


addition  of  the  damping  parameter  has  the  effect  of  reducing 
the  number  of  degrees  of  freedom,  the  trade-off  between 
resolution  and  standard-error  should  be  examined  in  order  to 
select  a  proper  choice  for  02.  For  this  study  we  chose 
02  *  0.1  (sec/%)  . 

Selection  of  the  initial  model  geometry  requires  experi¬ 
mentation  with  different  configurations  in  order  to  derive 
a  well-resolved  solution  with  small  standard  errors.  The 
blocksize  depends  on  the  station  distribution,  and  the  number 
of  observations.  Because  it  is  desirable  to  have  a  minimum 
of  ten  rays  sample  each  block,  we  chose  our  blocks  to  be  one 
degree  on  a  side.  The  inversion  should  be  carried  to  a 
depth  of  at  least  one  array  diameter  which  we  selected  to 
be  500  km.  As  suggested  by  Ellworth  (1977) ,  a  suitable  height 
to  width  ratio  is  1.5:1  which  we  found  to  be  satisfactory. 
Experimentation  with  displaced  layers  (in  an  effort  to  de¬ 
couple  blocks)  and  rotated  coordinate  systems  showed  little, 
if  any,  improvement  on  the  final  model. 

The  initial  model  used  for  the  ray  tracing  was  a  layered 
spherical  model  with  constant  velocities  in  each  layer. 
Velocity  differences  of  up  to  10%  do  not  substantially  alter 
the  sampling  of  the  blocks.  We  average  velocities  from  a 
Herrin  earth  model  to  derive  the  layer  velocities  (Table  2) . 

The  velocity  perturbations  and  the  number  of  observations 
for  blocks  in  each  layer  of  the  final  model  are  shown  in 
Figure  9.  In  these  figures,  negative  velocity  perturbations 
represent  relatively  lower  velocities.  We  only  included 
blocks  which  were  sampled  by  a  minimum  of  ten  rays  which 
mainly  eliminated  the  outermost  elements. 
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Modeling  of  the  crust  is  probably  one  of  the  most  crucial 
results  we  wish  to  obtain.  Unfortunately,  because  of  small 
incidence  angles,  most  rays  pass  through  the  crustal  block 
and  the  block  lying  directly  underneath.  This  makes  it 
difficult  to  decouple  the  effects  of  the  two  superposed  blocks. 
Resolution  could  be  improved  if  the  thickness  of  the  surface 
layer  were  increased.  However,  the  surface  blocks  would 
then  contain  a  large  and.  important  velocity  discontinuity 
(the  Moho)  and  information  regarding  crustal  structure  would 
be  lost. 

Problems  may  also  arise  because  of  the  uneven  station 
distribution.  If  rectangular  prisms  are  used  to  model  the 
crust,  a  number  of  stations  located  in  areas  with  potentially 
different  structures  may  be  included  within  the  same  element. 

To  remedy  this  problem  we  compute  velocity  perturbations  for 
each  individual  station  rather  than  for  a  given  block  config¬ 
uration.  It  was  found  that  using  this  approach  had  little 
effect  on  the  upper  mantle  layers  and  it  facilitated  a  means 
to  compute  a  three-dimensional  crustal  thickness  model  which 
is  described  below. 

Although  resolution  is  poor  and  standard  errors  are  large, 
the  velocity  perturbations  show  a  good  inverse  correlation 
with  the  average  station  residuals  (Figure  10) .  This 
further  supports  previous  conclusions  that  the  average  station 
residuals  are  mainly  dependent  on  crustal  structure,  and  may 
reflect  crustal  thicknesss  and  velocity  variations.  Iyer 
and  Healy  (1972)  successfully  used  a  general  relationship 
to  model  crustal  thickness  at  L&SA  from  the  distribution  of 


average 


v. 


e  station  residuals.-.-  ..They  a 
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geneous  crust  beneath  the  array  which  may  be  a  reasonable 
assumption  for  a  first  order  crustal  model  in  the  particular 
region  they  were  modeling.  However,  in  the  previous  section, 
we  concluded  that  the  two  orogenic  belts  in  the  northeastern 
U.S.  show  regional  differences  in  average  crustal  velocities 
and  thicknesses.  This  implies  that  we  must  account  for 
lateral  velocity  variations  in  order  to  use  the  average 
residuals  for  crustal  thickness.  Combination  of  the  velocity 
perturbations  with  the  average  station  residuals  allows  us 
to  construct  a  crustal  thickness  model.  Although  there  is 
a  trade-off  between  crustal  thickness  and  lateral  velocity 
differences  we  feel  that  the  method  is  justified  because 
independent  geophysical  data  suggests  the  dependence  of 
average  station  residuals  on  crustal  thickness  and  lateral 
velocity  variations. 

5V 

We  assumed  that  the  velocity  perturbations,  represent 
differences  f rota  an  average  crustal  velocity,  V  ,  It  should 
be  stressed  that  any  information  regarding  absolute  velocities 
in  the  three  dimensional  inversion  is  lost  and  the  technique 
we  are  using  is  therefore  an  approximation  and  represents  only 
a  first  order  crustal  model.  The  average  crustal  velocity, 

V^,  beneath  each  station  i  is  then 

vi  “  vc((^7,i/100)  +  vc  {9) 

and  the  crustal  thickness  is 

hi  *  Vi5i  +  ho  (10> 


where  h„  is  some  average  crustal  thickness  to  be  varied.  We 

chose  V  ■  6.5  km/s  and  h  *  38  km  and  calculated  the  crustal 
c  o 

thickness  and  average  crustal  velocity  maps  shown  in  Figure 
11.  Overall,  the  crust  beneath  the  Appalachian  province 
appears  to  be  3  km  thicker  than  that  of  the  Grenville. 

Average  crustal  velocities  differ  very  little  but  may  be 
slightly  lower  in  many  parts  of  the  Appalachians.  There  is 
a  region  of  greater  crustal  thickness  beneath  eastern  Vermont, 
New  Hampshire,  and  southern  Maine.  In  western  Vermont,  and 
western  Connecticut  there  axe  high  thickness  gradients  where 
the  crust  thins  to  the  west.  These  features  are  located  just 
to  the  east  of  the  Taconic  thrust  belt  and  in  western  Vermont, 
the  gradients  show  an  impressive  correlation  with  the 
serpentinite  belt.  Slight  crustal  thinning  is  observed  in 
eastern  Massachusetts,  and  it  is  not  clear  if  this  represents 
a  trend  which  may  continue  eastward  towards  the  outer  edge  of 
the  continental  margin.  Alternatively,  it  may  represent  the 
contact  between  the  central  orogenic  belt  and  the  eastern  belt 
described  earlier.  The  average  crustal  thickness  of  36  km 
in  northern  New  York  agrees  very  well  with  refraction  models 
of  Katz  (1955)  and  those  of  Aggarwal  in  Schnerk  et  al  (1976) . 

The  upper  mantle  (layer  2)  between  35  and  200  km  (Figure  9b) 
depth  is  characterized  by  a  northeast  trending  region  of  low 
velocities  beneath  Massachusetts ,  New  Hampshire  and  southern 
Maine.  The  trend  of  this  zone  has  an  interesting  parallel 
with  the  "grain”  of  the  Appalachian  surface  structure.  The 
velocities  appear  to  increase  to  the  west  and  northwest  be- 


neath  the  Precambrian  Grenville  Province.  This  suggests 
that  structural  differences  between  the  two  orogenic  belts 
extends  into  the  upper  mantle.  The  resolution  in  the  upper 
mantle  layer  increases  markedly  relative  to  that  of  the  crustal 
elements.  Because  the  block  configuration  in  layer  2  is 
better  suited  for  sampling  incident  waves,  each  element, 
particularly  the  central  elements,  observe  a  good  cross-fire 
of  rays.  Resolution  degrades  somewhat  towards  the  outer 
edges  of  the  model  because  the  number  of  observations  de¬ 
creases,  and  the  azimuths  cure  limited  to  a  small  sector. 

The  overall  resolution  and  standard  errors  are  best  for  this 
particular  layer. 

In  layer  3  (200-350  km)  the  region  of  lower  velocities  is 
shifted  west  and  northwest  and  is  found  beneath  Vermont  and 
central  Maine.  There  is  still  some  evidence  of  higher 
velocities  in  the  western-most  blocks  in  the  model. 

The  deepest  layer  (layer  4;  350-500  km)  shows  little 
resemblance  to  the  overlying  layers.  The  overall  trends 
appear  to  be  of  shorter  wavelength  and,  not  surprisingly, 
there  is  little,  if  any  distinction  between  the  two  geologic 
provinces  observed  at  the  surface. 


Discussion 


Three-dimensional  inversion  of  travel  time  data  illus¬ 
trates  that  structures  down  to  perhaps  200  km  and  greater 
can  be  correlated  with  surface  geology.  This  has  the  important 
implication  that  major  orogenic  belts  have  effects  that  reach 
well  into  the  lithosphere  which  are  stable  for  extended 
periods  of  time,  perhaps  as  long  as  1  billion  years. 

In  addition  to  the  regional  differences  existing  between 
the  two  structural  provinces  in  the  northeastern  U.S.,  we 
see  smaller  scale  features  that  are  related  to  specific 
segments  of  each  orogenic  belt  and  to  the  intervening  suture 
zone.  As  discussed  previously,  the  Taconite  belt,  the  east¬ 
ward-  lying  Precambrian  uplifts  (The  Berkshire  Highlands*  and 
Green  Mountains) ,  and  serpentinite  belt  are  Lower  Ordovician 
'  Structures  related  to  the  closure  of  an  inner-arc  basin  and 
a  continent-island  arc  collision.  This  region  is  not  only 
a  site  of  large  scale  geologic  deformation  but  it  is  also 
characterized  by  geophysical~ahomalies  such  as  apparent 
crustal  thinning,  high  crustal  velocities,  and  large  positive 
Bouguer  gravity  anomalies  (Figure  6) .  It  is  possible  that 
these  effects  are  related  to  deep-seated  thrusts  emplaced 
near  the  end  of  the  Taconic  orogeny  which  carried  oceanic 
crust  and  upper  mantle  material  to  higher  levels  within  the 
crust.  This  material  would  have  high  velocities  and  densities 
relative  to  the  surrounding  eugeoclinal  lithologies  and  would 
provide  a  model  consistent  with  the  observed  geophysical  and 
geological  anomalies.  There  are  at  least  two  objections  to 


this  model.  The  first  is  that  higher  crustal  velocities  are 
observed  just  to  the  west  of  the  Taconic  thrust  belt  at 
stations  PNY  and  WNY.  Secondly,  we  see  that  geophysical 
anomalies  are  not  observed  along  the  entire  Taconide  belt, 
although  the  geologic  structures  are  fairly  uniform.  This 
gap  in  the  zone  of  high  crustal  thickness  and  velocity 
gradients  east  of  and  parallel  to  the  thrust  belt  in  south¬ 
western  Vermont  and  western  Massachusetts  may  simply  be 
an  effect  of  poor  station  coverage  in  the  region.  As  pointed 
out  by  Fletcher  et  al.  (1978) ,  the  anomalous  residuals  beneath 
PNY  may  be  related  to  buried  mafic  or  ultramafic  bodies  whose 
presence  is  indicated  by  local  gravity  highs  (anomaly  22  of 
Simmons,  1964).  However,  it  is  readily  apparent  that  all 
of  the  stations  in  northern  New  York  record  very  early  arrivals 
from  the  northwest.  This  implies  the  existence  of  regional, 
deep-seated  effects  extending  to  the  northwest  beneath  the 
Adirondack  dome  and  is  supported  by  the  three-dimensional 
inversion  (Figure  9) .  Interestingly,  the  large  negative 
residuals  from  the  northwest  are  not  observed  at  stations 
further  south  which  suggests  that  the  zone  of  relatively 
high  upper  mantle  velocities  is  controlled  by  a  structure 
unique  to  the  Adirondacks.  It  has  been  postulated  that  the 
Adirondacks  and  associated  Grenville  basement  represent 
deeper  levels  of  an  ancient  "Tibetan  Plateau"  characterized 
by  crustal  thickening  and  shortening  behind  a  continental 
collision  zone  (Dewey  and  Burke,  1973;  Toksflz  and  Bird,  1977). 
It  is  not  clear  when  the  Adirondack  dome  was  elevated,  and 


the  cause  of  relatively  high  velocities  found  in  the  upper 
mantle  beneath  the  region  remains  problematical.  If  the 
Grenville  belt  does  represent  an  analog  of  a  Tibetan  Plateau, 
and  crustal  thicknesses  presently  average  36  km,  then  possibly 
another  15-35  km  of  crust  has  been  eroded  away  since  the 
late  Precambrian.  This  implies  tha t  the  presently  exposed 
surface  geologic  features  of  the  Adirondacks  represent 
lithologies  deformed,  metamorphosed,  and  intruded  at  great 
crustal  depths. 

Another  prominent  feature  in  the  crustal  model  is  the 
apparent  crustal  thickening  beneath  central  New  Hampshire  and 
southern  Maine.  Topographically,  these  are  the  regions  of 
highest  elevations  in  New  England  and  are  therefore  iso- 
statically  compensated  highland  areas.  The  Lower  Devonian 
Acadian  orogeny  had  a  severe  impact  on  the  structure  as 
evidenced  by  high  grade  metamorphism,  intrusion  of  granitic 
plutons,  and  large  scale  westward  recumbent  folding  of  many 
of  these  plutons  (Naylor,  1968).  Many  lines  of  evidence 
suggest  the  uplift  of  a  significant  mountain  belt.  Early 
Devonian  miogeoclinal  sediments  in  eastern  New  York  are 
overlain  by  a  thick  clastic  Middle  Devonian  sequence  of  the 
Catskill  Delta  -  a  thick  clastic  wedge  derived  from  the  up¬ 
lifted  landmass  in  central  New  England  (Rodgers,  1970).  The 
crustal  thinning  towards  central  and  northern  Maine  is 
correlated  with  the  drop-off  in  the  grade  of  metamorphism 
(Thompson,  1968) .  This  implies  that  the  continental  collision 
of  Acadian  orogeny  was  particularly  severe  in  southeastern  and 
central  New  England  and  we  would  expect  to  see  greater  crustal 
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thicknesses  in  these  regions. 

In  eastern  Massachusetts  we  again  see  an  apparent  crustal 
thinning  over  much  of  the  eastern  belt  described  in  an  earlier 
section.  These  contours  conform  very  well  with  the  position 
of  the  middle  Paleozoic  thrust  belt  and  the  abrupt  falloff 
in  metamorphic  isograds  from  sillimanite  to  chlorite 
(Thompson,  1968). 

The  low  velocity  anomalies  in  the  upper  mantle  (35-200  km) 
beneath  Massachusetts,  New  Hampshire,  and  southcentral  Maine 
strike  northeast  parallel  to  the  structural  grain  of  the 
Appalachian  orogenic  belt.  There  is  some  evidence  that  this 
feature  dips  to  the  northwest  to  depths  greater  than  200  km. 
Because  of  the  geometry  of  the  initial  model,  the  vertical 
extent  of  this  structure  is  only  constrained  to  lie  somewhere 
between  200  and  350  km.  It  is  interesting  to  note  the  spatial 
correlation  of  this  low  velocity  trend  with  the  Bronson  Hill  - 
Boundary  Mountains  anticlinorium.  As  discussed  previously, 
these  structures  probably  occupy  the  sites  of  a  complex  series 
of  island  arc  sequences  last  active  in  the  Early  Devonian 
time  prior  to  the  Acadian  orogeny.  If  it  is  true  that  sub- 
duction  occurred  east-southeast  of  the  Bronson  Hill  in  a 
west-northwest  direction,  we  may  expect  to  observe  some 
anomalous  feature  at  these  depths. 

Assuming  that  a  subducted  oceanic  slab  can  produce  such 
a  long-term  anomaly,  the  existence  of  velocities  that  are 
lower  than  the  surrounding  mantle  material  is  surprising, 
yet  not  necessarily  inconsistent  with  its  possible  thermal 
evolution.  Thermal  models  of  downgoing  oceanic  lithosphere 


(ToksdSz  et  al. ,  1973)  illustrate  that  in  the  early  stages 
o£  subduction,  the  slab  is  cold  and  dense  relative  to  its 
surroundings  to  a  depth  o£  at  least  600  km.  At  approximately 
100  m.y.  from  the  cessation  of  subduction,  the  slab  reaches 
thermal  equilibrium.  Thermal  models  suggest  that  beyond 
about  100  m.y.  the  subducted  oceanic  lithosphere  starts 
to  heat  up  because  of  its  higher  radioactivity  relative  to 
the  surroundings.  The  anomaly  is  not  necessarily  thermally- 
induced  and  other  explanations  related  to  the  compositional 
changes  of  the  descending  lithosphere  may  be  constructed. 

If  we  assume  the  velocity  anomalies  are  produced  by  temperature 
changes,  and  take  temperature  coefficient  of  velocity  for 
dunite,  a  1%  velocity  decrease  corresponds  to  a  temperature 
increase  of  about  100-200* C.  These  observations  are  con¬ 
sistent  with  preliminary  unpublished  heat  flow  measurements 
in  central  New  England  corrected  for  radioacitivity  which 
show  regions  of  highest  heat  flow  in  central  and  southeastern 
New  Hampshire  relative  to  northern  New  Ha  pshire  and  south¬ 
eastern  Massachusetts  (Jaupart,  pers.  comm.).  Geomagnetic 
data  from  the  northeastern  United  States  suggests  a  200° C 
temperature  anomaly  in  the  lower  crust  beneath  central  New 
Hampshire  and  southern  Maine  (Bailey  et  al . ,  1978). 

It  is  also  interesting  to  note  that  similar  trends  in 
teleseisaic  PP  residuals  are  observed  beneath  Newfoundland 
(Stewart,  1978) .  Newfoundland  is  a  northward  extension  of 
the  Appalachian  orogenic  belt  and  it  is  not  inconsistent  that 
the  crust  and  upper  mantle  structures  appear  to  be  similar 
to  our  observations  in  New  England. 


Figure  Captions 


Figure  la.  Generalized  geologic  map  of  northeastern  United 
States.  Modified  from  King  (1969). 

Figure  lb.  Major  tectonic  structures  in  northeastern  United- 
States.  Shaded  regions  correspond  to  sillimanite-grade 
metamorphic  zones. 

Figure  2.  Seismograph  stations  used  for  three-dimensional 
inversion  of  teleseismic  P  wave  residuals.  Labeled 
stations  correspond  to  those  referred  to  in  the  text 
or  shown  in  Figure  5. 

Figure  3.  Epicenters  of  earthquakes  and  explosions  used  as 
P  wave  sources. 

Figure  4.  Map  of  average  relative  travel  time  residuals  in 
seconds  for  seismograph  stations  in  northeastern  United 
States. 

Figure  5.  Station  focal  sphere  projection  plots.  Relative 

residuals  plotted  as  a  function  of  azimuth  to  the  source 
and  incidence  angle. 

Figure  6.  Bouguer  gravity  in  northeastern  United  States. 

Areas  with  positive  gravity  anomalies  are  shaded. 

Modified  from  Kane  et  al. ,  (1972). 

Figure  7.  Mean  relative  teleseismic  travel  time  residual  in 
seconds  relative  to  Bouguer  gravity  for  northeastern 
United  States  stations.  Open  circles  correspond  to 
stations  located  in  the  Grenville  Province  above  Pre- 
cambrian  basement.  Closed  circles  correspond  to  stations 
mainly  in  the  Appalachian  Province  located  above  Paleozoic 
basement . 


In  summary,  we  have  presented  geophysical  data  which 
suggests  that  structural  differences  between  two  juxtaposed 
orogenic  belts  exist  to  depths  of  200  tan  and  greater.  This 
indicates  that  major  collisional  episodes  have  very  long 
lasting  deep-seated  effects  as  well  as  large  areal  extents. 
Crustal  thickness  and  velocity  variations  are  correlative 
with  many  surficial  geologic  and  tectonic  features.  The 
Paleozoic  Appalachian  province  appears  to  have  a  slightly 
thicker  crust  with  lower  average  velocities  than  the  Pre- 
cambrian  Grenville  province.  The  higher  average  velocities 
associated  with  the  Grenville  province  are  evident  to  depths 
of  200  km  particularly  beneath  the  Adirondack  dome.  A 
relatively  low  velocity  anomaly  extending  to  depths  of  at 
least  200  km  and  dipping  to  the  northwest  shows  a  spatial 
correlation  with  the  Bronson  Bill  -  Boundary  Mountains 
Anticlinorium.  This  may  be  due  to  effects  of  post-Acadian 
radioactive  heating  of  subducted  oceanic  lithosphere. 
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Table  2 

Initial  Layered  Velocity  Model  For  Three-Dimensional  Inversion 
Layer  No.  Velocity (km/s)  Thickness (km)  Block  Length 
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Normal  fault,  hachures  on  downthrown  side 
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ABSTRACT 


The  attenuation  of  seismic  waves  in  upper  crustal 
rocks  is  studied  from  both  experimental  data  and  theoretical 
models  of  specific  attenuation  mechanisms.  The  attenuation, 
or  anelasticity  may  be  used  in  conjunction  with  the 
compressional  (P)  and  shear  (S)  wave  velocities  to  infer 
the  microstructure  of  the  rock  and  the  physical  conditions 
encountered  in  the  subsurface  environment.  Generally, 
the  attenuation  varies  much  more  than  velocity  with 
accompanying  changes  in  the  physical  state  making  this 
measurement  a  particularly  valuable  tool.  Unfortunately, 
systematic  studies  of  the  behavior  of  attenuation  with 
parameters  such  as  pressure  and  saturation  conditions  are 
deficient.  Nor  has  the  development  of  new  experimental 
techniques  and  theoretical  models  kept  pace  with  the 
methods  for  elastic  properties. 

It  is  shown  in  this  thesis  that  previous  experimental 
data  suggest  one  or  more  of  several  proposed  mechanisms 
may  be  contributing  to  the  overall  attenuation  observed 
in  crustal  rocks.  These  include  crack  and  grain  boundary 
sliding  friction;  fluid  associated  mechanisms  such  as 
crack  lubrication,  fluid  flow  through  porous  and 
permeable  rock,  squirting  flow  from  thin  cracks  to  larger 
pores,  and  viscous  shear  relaxation  in  the  void  spaces; 
losses  associated  with  gas  bubble  resonances;  and  a  host 
of  geometrical  effects  with  scattering  from  pores  in  the 
rock  being  most  important  from  an  experimental  point  of 
view.  Mathematical  models  of  all  these  mechanisms  in 
terms  of  their  frequency  and  pressure  dependences  are 
developed. 

The  attenuation  or  Q  values  (inverse  proportional  to 


the  energy  loss)  for  P  and  S  waves  are  determined  for  a 
suite  of  upper  crustal  rocks  as  functions  of  pressure 
and  saturation  using  an  ultrasonic  pulse  transmission 
technique  with  the  goal  of  interpreting  these  results  in 
terms  of  the  theoretical  models  of  the  mechanisms.  In 
general ,  it  is  found  that  Q  for  both  P  and  S  waves 
increases  (attenuation  decreases)  as  a  function  of  pressure 
for  dry  and  Water  saturated  rocks.  Q  values  for  saturated 
and  partially  saturated  rocks  are  less  than  for  dry  rocks 
and  while  not  conclusive,  the  experiments  suggest  that 
Qp  for  low  gas  saturations  is  slightly  lower  compared  to 
the  fully  water  saturated  case.  The  major  controlling 
factor  in  determining  the  differences  in  Q  for  individual 
rocks  is  the  nature  and  concentration  of  thin  cracks  - 
the  exceptions  studied  being  shales  and  tuffs.  Rocks 
with  a  higher  density  of  fine  cracks  exhibit  lower  Q  values 
and  the  rate  of  change  of  Q  with  pressure  may  be  correlated 
with  the  crack  width  distribution  inferred  from  static 
stress-strain  measurements. 

Modeling  of  the  Q  data  obtained  for  dry  and  water 
saturated  samples  of  Berea  sandstone  imply  that  friction 
is  the  dominant  mechanism  for  attenuation.  Increasing 
pressure  decreases  the  number  of  sliding  surfaces  and  thus 
Q  increases.  The  major  effect  of  fluid  is  to  lubricate 
these  surfaces  thus  facilitating  sliding.  Attenuation  due 
to  fluid  flow  plays  a  secondary  role.  At  ultrasonic 
frequencies  and  moderate  pressures,  flow  is,  however, 
important  in  porous  and  permeable  rock  like  the  sandstone. 
Scattering  is  dominant  only  at  very  high  frequency  or,  as 
with  the  tuffs,  when  the  grain  or. poire  sizes  are  large 
compared  to  the  wavelength.  Attenuation  in  shales  may  be 
due  to  relative  shear  motion  between  the  rock  matrix  and 
clay  or  kerogen.  Bubble  associated  losses  in  partially 
saturated  rocks  may  be  important  only  in  limited  cases. 

The  phenomenon  of  frictional  attenuation  is  further 
examined  using  a  resonance  bar  experimental  technique. 

It  is  found  that  Q  in  dry  rocks  is  dependent  on  amplitude, 
constant  at  low  strains  but  decreasing  rapidly  at  strains 
greater  than  about  10-s,  explaining  differences  in  Q 
observed  between  resonance  and  ultrasonic  methods.  The 
amplitude  dependence  is  also  a  function  of  the  types  of 
cracks  existing  in  the  rock.  This  is  investigated  for 
samples  that  have  been  thermally  cycled  in  order  to  induce 
cracking.  It  is  found  that  rocks  cycled  at  low  heating 
rates  show  an  increase  in  Q  up  to  a  cycled  temperature  of 
about  400°C  which  is  accompanied  by  a  migration  of  the 
amplitude  transition  to  lower  strains.  These  data  are 
interpreted  in  terms  of  the  widening  of  thin  cracks  thus, 
reducing  the  contact  area  of  frictional  sliding  -  both 


increasing  Q  and  allowing  the  surfaces  to  slide  at  lower 
amplitudes.  Some  of  the  effect  observed  may  also  be  due 
to  outgassing  of  volatiles  from  the  crack  surfaces. 

The  results  obtained  in  this  thesis  may  be  used  to 
evaluate  the  difficulties  involved  in  extrapolating 
laboratory  data  to  insitu  measurements  of  attenuation. 
Theoretical  models  such  as  those  used  to  describe  the 
anelastic  behavior  of  the  Berea  sandstone,  provide  the 
best  method  at  the  present  time.  However,  field  measure 
ments  must  be  corrected  for  other  factors  that  affect 
seismic  wave  amplitude.  Furthermore,  the  contributions 
of  mechanisms  that  may  be  important  at  low  frequencies 
and  not  easily  established  from  ultrasonic  data  plus  the 
amplitude  dependence  of  attenuation  must  be  considered 
before  any  meaningful  comparisons  can  be  made. 
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ABSTRACT 


The  role  of  crack  and  grain  boundary  contacts  in 

determining  seismic  wave  attenuation  in  rock  is  investigated 

by  examining  Q  as  a  function  of  thermal  cycling  (cracking) 

and  wave  strain  amplitude.  Q  values  are  obtained  using  a 

longitudinal  resonant  bar  technique  in  the  10-20  kHz  range 
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for  maximum  strain  amplitudes  varying  from  roughly  10  to 

10  The  samples  studied  include  the  Berea  and  Navajo 

sandstones,  Plexiglas,  westerly  granite,  Solenhofen 

limestone  and  Frederick  diabase,  the  latter  two  relatively  crack 

free  in  their  virgin  state.  Measurements  were  made  at  room 

temperature  and  pressure  in  air. 

Q  values  for  both  sandstones  are  constant  at  low  strains 

C<10  ).  but  decrease  rapidly  with  amplitude  at  higher  strains. 
There  is  no  hysteresis  of  Q  with  amplitude.  Q  values  for 

Plexiglas  show  no  indication  of  amplitude  dependent  behavior. 

The  granite,  limestone  and  diabase  are  thermally  cycled 

at  both  fast  and  slow  heating  rates  in  order  to  induce 

cracking.  Samples  slowly  cycled  to  400°C  show  a  marked 

increase  in  Q  that  cannot  be  entirely  explained  by  outgassing 

of  volatiles.  Cycling  widens  thin  cracks  and  grain  boundaries. 


reducing  contact  areas.  Rapidly  heated  samples  exhibit 
completely  different  behavior.  Amplitude  dependence  of  Q 
is  found  to  be  coupled  to  the  effects  of  thermal  cycling. 
The  transition  from  low  to  high  strain  amplitude  behavior 


occurs  at  lower  strains  for  thermally  cycled  rocks. 


INTRODUCTION 


The  presence  of  fine,  thin  cracks  and  grain  boundary 
contacts  has  been  shown  to  have  an  important  role  in 
determining  many  of  the  elastic  and  anelastic  properties  of 
rock.  In  particular,  they  control  the  velocity  and  attenuation 
behavior  at  moderate  confining  pressures  (Walsh,  1965  and  1966; 
Toksfiz  et  al.,  1976;  Johnston  and  Toksflz,  1979).  Increases 
in  both  velocity  and  Q  with  pressure  at  all  frequencies  and 
under  a  variety  of  saturation  and  strain  amplitude  conditions 
may  be  attributed  to  the  closing  of  these  cracks  and  grain 
boundary  contacts.  The  pressure  dependence  of  Q,  however, 
provides  only  a  general  indication  as  to  the  importance  of 
cracks.  As  will  be  seen,  attenuation  does  not  depend  simply 
on  parameters  like  crack  density  but  rather  in  more  complicated 
ways  on  distributions,  shapes,  sizes,  and  contacts. 

The  variation  of  Q  with  wave  strain  amplitude  also  proves 

to  be  an  important  diagnostic  for  the  understanding  of  contact 

related  attenuation.  It  is  almost  invariably  assumed  in  the 

literature  that  attenuation  is  independent  of  strain  amplitude. 

While  this  may  indeed  be  true  for  strains  associated  with 

teleseismic  waves  (£<10“®)  numerous  experimental  studies  have 

shown  amplitude  dependent  behavior  at  higher  strains  (£>10_<*). 

Gordon  and  Davis  (1968)  examined  Q  as  a  function  of  amplitude 
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ranging  from  10  to  about  10  using  a  resonant  bar  at  90  kHz. 
The  samples  studied  included  granite,  basalt,  dunite,  quartzite 
and  others.  Q  was  observed  to  be  constant  at  low  strains. 


decreasing  sharply  for  strain  amplitudes  greater  than  about 
10  Other  studies  by  Peselnick  and  Outerbridge  (1961)  on 
Solenhofen  limestone  under  torsional  vibrations  at  1  Hz  and 
Gordon  and  Rader  (1971)  for  the  Chester  granite  and  Winkler 
and  Nur  (1978)  for  the  Massilon  sandstone  produce  similar  results. 

As  will  be  shown,  the  observation  of  amplitude  dependence 
at  strains  of  about  10-^  and  greater  may  also  explain  many  of 
the  discrepancies  reported  between  Q  measurements  obtained  by 
dynamic  resonance  and  ultrasonic  pulse  techniques  although  a 
frequency  dependent  attenuation  mechanism  cannot  be  ruled  out. 

To  investigate  the  nature  of  amplitude  dependent  Q,  we 
carried  out  two  sets  of  experiments.  In  the  first,  the  amplitude 
dependence  of  Q  was  investigated  for  two  sandstones  and  Plexiglas. 
In  the  second,  the  behavior  of  attenuation  was  examined  as  a 
function  of  cracks  induced  by  temperature  cycling  of  the  rock. 

For  this  case,  the  effects  of  both  increasing  crack  density 

and  outgassing  of  volatiles  from  crack  surfaces  must  be  considered. 


Furthermore,  it  was  observed  that  amplitude  dependence  is  coupled 
to  the  effects  of  thermal  cycling. 


EXPERIMENTAL  RESULTS 

Q  values  for  all  the  experiments  reported  in  this  paper 
were  determined  using  a  longitudinal  resonant  bar  method, 
driven  as  a  composite  oscillator  (Quimby,  1925)  by  high  Q  PZT-4 
transducers  cemented  to  the  sample  ends.  The  samples  were 
about  15  cm  long  by  0.635  cm  in  diameter  producing'  fundamental 
mode  frequencies  in  the  10-20  kHz  range.  Several  runs  were 
done  with  a  sample  diameter  of  1.27  cm.  Bar  velocity  and  Q 


are  found  from  the  resonance  frequency  and  half  power 
bandwidth.  These  values  are  corrected  for  the  loading  effects 
and  mechanical  Q  of  the  transducers  (Bozorth  et  al.,  1951; 

Marx,  1951) . 

All  of  the  measurements  were  made  under  room  conditions 

with  unjacketed  samples.  Maximum  strain  amplitude  at  resonance 
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was  varied  from  roughly  10  to  10  by  increasing  the  peak 
voltage  to  the  driving  transducer. 

The  samples  studied  include  the  Berea  and  Navajo  sandstones, 
Plexiglas,  Westerly  granite,  Frederick  diabase,  and  Solenhofen 
limestone.  The  latter  two  were  chosen  for  their  low  crack 
porosity  in  the  virgin  state.  Descriptions  of  the  sandstones 
may  be  found  in  Johnston  and  Toksdz  (1979),. the  granite  and 
diabase  in  Brace  (1965)  and  the  limestone  in  Birch  (1960) . 

As  mentioned  in  the  introduction,  attenuation  values  for 
dry  rocks  obtained  by  dynamic  resonace  and  ultrasonic  techniques 
differ,  the  former  yielding  higher  Q's.  Since  these  measurements 
are  generally  made  at  different  frequencies,  one  might  possibly 
hypothesize  a  frequency  dependent  mechanism.  However,  no  data 
support  this.  Q  dependent  of  strain  amplitude  though,,  may 
provide  an  explanation.  Most  resonance  experiments  are  carried 
out  at  low  strains  relative  to  those  associated  with  ultrasonic 
pulses.  If  Q  for  low  strains  is  higher  than  for  high  strains, 
the  apparent  discrepancy  in  the  data  may  be  resolved.  This 
behavior  is  indeed  observed  for  the  dry  Berea  sandstone  as 
shown  in  Figure  1  where  Q  and  bar  velocity  are  plotted  as  a 
function  of  maximum  strain  amplitude.  Also  shown  for  comparison 


is  Q  for  amorphous  Plexiglas  over  the  same  amplitude  range. 
Behaving  similarly  to  data  reported  by  Winkler  and  Nur  (1978) 
for  the  Massilon  sandstone,  Q  for  the  Berea  rapidly  decreases 
from  a  value  of  75  to  a  value  of  30  at  a  strain  amplitude  of 
about  10"5.  This  corresponds  to  a  60%  decrease  in  Q.  Velocity 
decreases  by  only  1%.  Also  note  that  both  velocity  and  Q 
appear  to  level  off  at  higher  amplitudes. 

Comparison  of  the  increasing  amplitude  curve  versus 
decreasing  amplitude  shows  nearly  complete  recovery  of  the  initial 
Q  value.  As  will  be  seen,  this  lack  of  hysteresis  is  observed  for 
all  the  rocks  studied  and  implies  that  little  or  no  damage  to  the 
sliding  surfaces  (assuming  that  attenuation  is  due  to  frictional 
sliding)  occurs  at  the  higher  strains. 

For  the  same  range  of  strain  amplitude,  Q  of  the  Plexiglas 
remains  relatively  constant  and  if  anything,  increases  slightly. 
This  fact  strongly  suggests  that  the  presence  of  grain 
boundaries  and  cracks  in  the  rock  is  the  controlling  factor  in 
producing  amplitude  dependent  behavior. 

In  Figure  2,  Q  is  plotted  as  a  function  of  strain  amplitude 
for  the  Navajo  sandstone.  The  same  behavior  as  seen  for  the 
Berea  -  constant  Q  at  low  amplitudes  with  a  sharp  decrease  at 
a  strain  amplitude  of  10  ®  is  observed,  although  the  transition 
from  low  to  high  amplitude  behavior  does  seem  to  occur  at 
slightly  lower  amplitudes  in  the  Navajo. 

Does  the  amplitude  dependence  explain  the  difference  between 
reported  resonance  and  ultrasonic  measurements?  For  the  most  part 
-  yes.  Qp  values  at  low  pressure  reported  by  Toksftz  et  al.  (1979) 
and  Johnston  and  ToksiJz  (1979)  using  ultrasonic 


pulses  are  compared  to  the  results  just  described  for  the 
dry  Berea  and  Navajo  sandstones  in  Table  1.  Here,  low  and 
high  amplitude  resonant  bar  Q  values  are  listed  along  with  the 
ultrasonic  P  wave  values.  Clearly  the  high  amplitude  results 
from  the  resonance  technique  are  comparable  to  the  ultra¬ 
sonic  data.  The  data,  however,  are  limited  and  it  would  be 
premature  to  rule  out  other  mechanisms-  even  frequency 
dependent  dry  friction  -  as  supplements  to  amplitude 
dependence  in  interpreting  the  discrepancy  between  the  two 
experimental  methods.  It  is  inescapable,  however,  that 
amplitude  dependence  must  be  considered  if  any  comparison 
is  to  be  made  between  not  only  different  experimental  tech¬ 
niques  but  also  between  laboratory  and  field  data. 

-he  preceding  daca  stronclv  sucjqcet  that  cracr.e  in  the 
rock  are  important  in  determining  the  amplitude  dependent 
behavior  of  Q.  Thermal  cycling  presents  a  unique  opportunity 
to  form  cracks  in  rock.  Thus,  if  information  on  the  induced 
crack  porosity  is  available,  the  change  in  Q  might  be 
quantified  in  terms  of  crack  density.  Intuitively,  one 
would  expect  that  Q  would  decrease  with  increasing  thermal 
cycling.  As  will  be  seen,  this  is  not  always  the  case.  In 
fact,  one  of  the  more  interesting  results  is  what  can  be 


said  about  the  nature  of  thermal  cracking  mechanisms. 

Samples  of  Frederick  diabase,  Solenhofen  limestone, 
and  Westerly  granite,  0.635  cm  in  diameter,  were  thermally 
cycled  up  to  a  maximum  temperature  of  600 °C  at  room  pressure 
in  order  to  induce  cracking.  Second  samples  of  each  rock 
were  cycled  in  separate. runs  in  order  to  check  reproducibility. 
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The  heating  and  cooling  rates  were  kept  under  28C/minute. 

The  samples  were  exposed  to  the  maximum  temperature  for  at 
least  one  hour.  Data  from  Johnson  et  al.  (1978)  suggest 
that  this  time  is  required  for  the  decay  of  acoustic  emissions 
(and  thus  cracking)  in  a  rock  at  constant  temperature. 

Other  samples  of  each  rock  were  slowly  cycled  to 
maximum  temperatures  of  400°  and  800 °C.  These  were  used 
for  static  stress-strain  measurements  to  obtain  a  handle 
on  the  induced  crack  density.  Finally,  samples  of  Solenhofen 
limestone  and  Frederick  diabase,  1.27  cm  in  diamater,  were 
subjected  to  thermal  cycling  at  high  heating  rates  of  about 
20°C/minute.  The  samples  were  maintained  at  the  maximum 
temperature  for  at  least  an  hour  to  equalize  outgassing 
effects  compared  to  the  slow  heating  rate  experiments. 

Westerly  granite  did  not  remain  coherent  during  these  runs 
and  was  thus  not  used. 

In  both  procedures,  after  cooling  the  samples  were  left 
in  laboratory  air  for  several  days.  Again  note  that  all 
attenuation  measurements  were  made  at  room  temperature  and 
pressure  in  air.  The  relative  humidity  in  the  laboratory 
varied  during  the  course  of  the  experiments  from  45  to  75%. 
This  appeared  to  have  had  no  bearing  on  the  Q  measurements. 

In  fact,  several  samples  were  saturated  with  water  and  then 
allowed  to  dry  with  little  change  in  the  results. 

Thermally  induced  crack  porosity  may  be  characterized 
using  static  stress-strain  analysis.  Linear  strain  as  a 
function  of  hydrostatic  confining  pressure  is  shown  for 
virgin  and  thermally  cycled  samples  of  diabase,  limestone, 
and  granite  in  Figures  3  through  7.  For  Westerly  granite. 


data  is  shown  for  two  orientations:  the  stiff  (parallel) 
and  soft  (perpendicular)  directions.  The  induced  crack 
porosities  for  cycling  to  400  and  800°  are  found  using  the 
method  described  by  Walsh  (1965)  and  the  results  are  listed 
in  Table  2.  The  increase  in  porosity  is  generally  propor¬ 
tionally  greater  for  the  samples  cycled  to  800°  than  those 
heated  to  400®.  The  shapes  of  the  thermally  induced  cracks 
seem  also  to  depend  on  the  temperature  achieved  during 
cycling.  In  the  diabase,  for  the  400° curve,  the  cracks  appear 
to  be  of  fairly  high  aspect  ratio.  The  closure  pressure  is 
relatively  high  (>1000  bars)  and  there  is  little  curvature 
at  low  pressures.  The  case  for  higher  aspect  ratio  cracks  in 
the  400° cycled  Westerly  samples  is  not  nearly  as  clear  as 
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aspect  ratio  cracks.  However,  evidence  to  be  presented 
later  suggests  relatively  higher  concentrations  of  wide  cracks 
compared  to  thin  ones  produced  by  thermal  cycling  to  temper¬ 
atures  no  greater  than  400® C.  In  both  the  diabase  and  granite, 
heat  treatment  to  800®  not  only  produces  a  dramatic  increase 
in  crack  porosity  but  also,  these  cracks  seem  to  be  of  lower 
aspect  ratio.  In  both  cases,  the  curvature  of  the  stress- 
strain  curve  at  low  pressures  is  strikingly  high  compared  to 
the  other  runs.  For  the  Solenhofen  limestone,  however,  a 
rapid  change  in  crack  porosity  is  not  observed  after  heating 
to  800°  nor  are  the  induced  cracks  of  particularly  low 
aspect  ratio. 

Q  values  at  low  strain  amplitude  and  bar  velocity  as 
functions  of  maximum  temperature  obtained  in  slow  heating 


and  cooling  cycles  for  the  diabase  and  limestone  are  shown 
in  Figures  8  and  9.  Q  and  V  for  the  two  core  directions  of 
Westerly  are  shown  in  Figure  10.  For  all  three  rocks,  two 
separate  experiments  were  run  to  check  reproducibility  (with 
the  exception  of  the  "stiff"  core  of  Westerly  granite) . 

These  are  denoted  on  the  figures  by  different  symbols.  While 
the  experiments  were  run  under  different  humidity  conditions, 
there  is  quite  satisfactory  agreement  between  the  two  runs 
for  all  three  rocks. 

The  remarkable  aspect  of  the  data  is  the  dramatic  in¬ 
crease  in  Q  for  all  the  samples  at  maximum  temperatures  of 
400°  and  less.  Simultaneously,  velocity  decreases  -  quite 
substantially  for  the  Westerly  granite  which  contains  more 
vrsvkir  initially.  Tht*  changes  in  velocity-  itr-an  I or  Li:e 
granite  are  nearly  identical  to  those  reported  by  Johnson 
et  al  (1978)  in  a  similar  experiment  on  Westerly  (Q  values 
were  not  obtained) .  Anisotropy  in  Westerly  is  apparent, 
amounting  to  10%  for  velocity  and  25%  for  Q.  The  velocity 
changes  in  the  diabase  and  limestone  are  low  compared  to  the 
granite  for  temperatures  below  400°  followed  by  a  rapid 
change  above  that  point.  This  is  reflected  by  the  changes 
in  crack  porosity  upon  thermal  cycling  listed  in  Table  2. 

It  is  clear  then,  from  all  the  data,  that  Q  increases  up 
to  a  point,  even  though  crack  porosity  is  also  increasing. 

This  change  in  Q  is  also  apparently  independent  of  orientation 
as  shown  for  the  Westerly  granite  in  Figure  10. 

Between  400  and  500°,  however,  the  trend  in  Q  reverses 
and  in  two  of  the  three. cases,  Q  decreases  substantially. 


These  data  are  very  similar  to  those  reported  by  Todd  et  al. 
(1972)  for  Westerly  granite  and  Fairfax  diabase.  The  high 
temperature  Q  values  of  the  limestone  do  not  follow  the  form 
of  the  other  two  rocks,  showing  only  a  slight  decrease. 

This  difference  might  have  been  expected  from  the  static 
measurements  where  the  Solenhofen  was  the  only  sample  not 
exhibiting  a  large  increase  in  low  aspect  ratio  cracks  at 
800° .  It  should  be  noted  that  the  decrease  in  Q  above  400° 
observed  for  the  diabase  and  granite  is  not  initially  due 
to  the  a- 6  transition  in  quartz  which  does  not  occur  until 
573°  C. 


Q  and  bar  velocity  for  the  1.27  cm  diamater  cores  of 
rapidly  heated  limestone  and  diabase  are  shown  in  Figures  11 
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ature  for  the  limestone  is  clearly  opposite  that  observed 
for  low  heating  rates.  For  the  diabase,  there  is  an  initial 
decrease  in  Q  at  200°  followed  by  a  gradual  rise  to  a 
maximum  at  500°  although  the  Q  values  at  this  point  are  far 
below  those  obtained  by  slow  heating  even  though  the  sample 
here  was  kept  at  temperature  for  at  least  one  hour  as  before. 
A  sharp  decrease  in  Q  is  noted  between  500  and  600*  C  probably 
due  to  the  phase  transition  in  quartz.  In  cycling  to  800° C, 
a  marked  change  in  the  sample  appearance  was  noted.  This  was 
accompanied  by  a  further  decrease  in  Q  and  a  very  rapid 
decrease  in  velocity. 

These  results  imply  that  a  different  crack  producing 


mechanism  is  operative  for  high  heating  rates.  Compared  to 
the  slow  rate,  one  might  expect  that  finer  (lower  aspect 


ratio)  cracks  are  generated.  The  data  are  more  conclusive 
for  the  limestone  where  not  only  Q  decreases  continuously 
but  the  decrease  in  velocity  is  larger  than  for  the  slow 


heating  rate  case.  For  a  given  volume  concentration, 
lower  aspect  ratio  cracks  produce  larger  changes  in  the 
elastic  moduli  compared  to  high  aspect  ratio  cracks  (Toksoz 
et  al. ,  1976).  Apparently  in  the  diabase,  permanent  out- 
gassing  of  volatiles  has  a  compensating  effect  on  the 
generation  of  thin  cracks  at  temperatures  between  300  and 
500° C.  The  minor  contribution  of  this  same  effect  may  be 
observed  in  the  limestone  where  Q  levels  off  between  400 
and  600° C. 

Amplitude  dependent  behavior  is  also  coupled  to  the 
nature  of  thermally  induced  cracks.  O  as  a  function  of 
maximum  strain  amplitude  is  shown  for  several  temperatures 
each  for  the  three  slowly  cycled  rocks  in  Figures  13  through 
15.  Amplitude  dependence  not  readily  apparent  in  the 
virgin  states  of  the  diabase  and  limestone  appears  for 
thermally  cycled  samples.  The  important  feature  to  note 
from  the  data  is  that  while  Q  increases  as  a  function  of 
temperatures  for  low  strains,  the  transition  from  low  to 
high  strain  amplitude  behavior  migrates  to  lower  amplitudes. 
It  also  appears  that  the  Q  values  for  high  strains  tend  to 
converge.  Once  the  maximum  temperature  exceeds  400  *C  and 
the  Q's  at  low  strains  decrease,  the  amplitude  transition 
moves  back  to  higher  strains.  The  exception  to  this  is  the 
Solenhofen  limestone.  Again  it  is  apparent  in  this  case 
that  the  nature  of  the  cracks  obtained  at  high  temperature 
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are  different  than  in  the  granite  or  diabase. 

Amplitude  dependent  behavior  was  also  observed  for  the 
rapidly  heated  samples  although  for  a  given  maximum  obtained 
temperature,  the  transition  seems  to  occur  at  a  higher 
strain  than  for  the  slowly  heated  samples.  This  is  shown 
in  Figure  16  for  the  Frederick  diabase  where  Q  as  a  function 
of  amplitude  is  plotted  for  the  virgin  rock  and  the  two  heat 
treated  runs.  The  amplitude  dependence  for  the  20° C/minute 
treated  sample  is  not  as  clearly  defined  as  for  the  other 
run.  Either  the  transition  does  occur  at  higher  strains  or 
the  amplitude  dependence  is  smeared  out  or  reduced  somewhat. 

The  behavior  of  amplitude  cycling  on  Q  is  nearly  iden¬ 
tical  for  all  the  rocks  including  the  sandstones  shown  in 
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little  or  no  damage  to  contact  surfaces  at  high  strain  amplitudes 

DISCUSSION  AND  CONCLUSIONS 

The  data  presented  in  the  preceding  section  may  be 
interpreted  in  terms  of  the  mechanisms  for  volatile  out- 
gassing  from  crack  surfaces,  thermally  induced  cracking  and 
amplitude  dependent  attenuation.  The  theoretical  basis  for 
the  latter  is  not  well  developed,  the  difficulty  being 
that  the  specifics  of  such  amplitude  behavior  is  very  model 
dependent.  Thus,  the  discussion  of  this  aspect  must  be 
qualitative.  The  increase  in  Q  observed  in  slow  thermal 
cycling  may  be  due  to  one  or  both  of  the  following  mechanisms. 

1)  Outgassing  of  volatiles  tied  to  mineral  grains  and 
crack  contact  surfaces.  If  attenuation  is  due  to  some  kind 


of  friction  mechanism,  this  has  the  effect  of  increasing  the 
friction  coefficient  thus  decreasing  the  attenuation  (increasing  Q) 
Tittmann  et  al.  (1974  and  1976) ,  in  their  exhaustive  survey  of 
the  effects  of  volatiles  on  the  acoustic  properties  of  lunar 
rocks,  have  shown  that  very  high  Q  values  (>3000)  can  be 
obtained  in  the  laboratory  by  repeated  thermal  cycling  and  the 
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application  of  hard  vacuum  (<10  torr) .  However,  outgassed 
samples  exposed  for  30  minutes  to  saturated  vapors  of  several 
volatiles,  similar  to  water  in  their  physical  properties,  at 
ambient  pressure  resulted  in  40-60%  decreases  in  Q  (Tittmann 
et  al. ,  1976) . 

Clearly,  the  major  question  at  this  point  is  whether 
exposure  of  the  samples  to  air  after  cycling  has  resulted  in 
substantial  reabsorption  of  water  vapor.  It  is  probable  that 
for  the  most  part,  it  has.  If  the  effect  of  outgassing  is  to 
lock  contact  surfaces,  then  the  migration  of  the  amplitude 
dependence  to  lower  amplitudes  would  not  be  observed.  Furthermore, 
if  outgassing  had  a  major  effect  on  the  Q,  then  this  should  be 
seen  in  the  rapidly  heated  samples  also.  Rather,  Q  decreased  or 
remained  relatively  constant  with  maximum  temperature.  The 
slight  rise  in  Q  observed  for  the  diabase  may,  however,  indicate 
that  part  of  the  overall  effect  is,  in  fact,  due  to  volatile 
loss. 


thermal  cyclinc 


The  overall  crack  density  must  increase 


during  cycling,  resulting  in  decreases  in  measured  bar 
velocities.  Thermal  cracking  for  low  temperature  gradients, 
though,  is  due  to  thermal  expansion  coefficient  mismatches 
across  grains  (Myklestad,  1942;  Edwards,  1951;  Eshelby,  1957). 
Stresses  induced  by  this  differential  thermal  expansion  are 
concentrated  between  grains,  the  effect  being  to  widen 
pre-existing  crack  contacts  and  grain  boundaries.  Thus, 
contact  surfaces  are  reduced  and  attenuation  decreases.  At 
higher  temperatures  and  high  heating  rates  new  cracks  are 
formed  and  propagated. 

There  is  some  evidence  to  support  the  contention  that 
expansion  and  contraction  of  grains  at  thermal  equilibrium 
has  the  etrect  or  opening  cracxs.  uirrerenpiax  scram  auaxysi 
of  several  rock  types  including  Westerly  granite  and  Frederick 
diabase  (Cooper  and  Simmons,  1977)  show  that  the  peak  in 
crack  porosity  as  a  function  of  closure  pressure  shifts  to 
higher  closure  pressure  and  thus  higher  aspect  ratio  cracks 
as  a  result  of  thermal  cycling.  Wong  and  Brace  (1978)  have 
shown  that  the  pressure  required  to  produce  reversible  thermal 
expansion  also  increases  after  cycling.  These  effects  have 
been  observed  qualitatively  in  Figures  3  through  6  showing 
linear  strain  as  a  function  of  pressure.  Examples  of  crack 
widening  due  to  slow  thermal  cycling  have  also  been  shown 
using  scanning  electron  microscope  photographs  (Sprunt  and 
Brace,  1974).  The  effect  of  thermal  cycling  to  400°C  in 


95 


Chelmsford  granite  is  to  widen  pre-existing  cracks  and  to 


form  new,  fairly  wide,  cracks  at  the  grain  boundaries.  Sprunt 
and  Brace  also  show  that  upon  thermal  cycling  there  is  a 
decrease  in  the  number  of  low  aspect  ratio  cracks  ( 1 0“ 4 ) 

while  a  substantial  increase  in  the  number  of  cracks  with 

_  2 

aspect  ratios  of  10  in  Westerly  granite. 

Thus,  wile  the  effect  of  cycling  a  rock  to  400°C  is  to 
increase  the  overall  crack  porosity,  these  seem  to  be  of  high 
aspect  ratio.  This  is  most  evident  for  the  diabase  and 
limestone,  more  or  less  crack  free  in  their  virgin  states. 

Further  heating  produces  a  dramatic  change  in  the  crack  porosity 
for  the  diabase  and  granite.  It  is  possible  that  above  400°C, 
the  rocks  are  unable  to  acomodate  any  more  thermal  expansion 
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I. 

Thus  we  see  a  sharp  increase  in  strain  for  low  confining  pressures. 
The  change  for  the  Solenhofen  at  800°C  (Figure  4)  is  not  nearly 
as  dramatic  and  this  is  reflected  in  the  Q  values  which  did 
not  suffer  a  dramatic  decrease  for  maximum  temperatures  above 
400°C.  The  limestone  probably  deforms  plastically  at  these 
temperatures  and  thus  fine  crack  growth  is  inhibited. 

For  high  heating  rates ,  local  thermal  gradients  in  the 
rock  may  cause  crack  growth  (Goodier  and  Florence,  1964; 

Todd,  1973;  Wong  and  Brace,  1978).  Calculated  maximum 
gradients  obrained  in  the  high  heating  rate  experiments 
exceed  a  critical  value  observed  by  Todd  (1973)  for  gradient 


induced  cracking.  Thermal  stresses  in  this  case  exist  over 
a  volume  larger  than  the  grains  and  thus  cracks  may  even 
propagate  through  grains,  increasing  contact  area. 


Strain  amplitude  effects  may  be  incorporated  into 
this  model  via  several  mechanisms.  Mavko  (1978)  has 
hypothesized  that  observed  amplitude  dependent  attenuation 
behavior  is  due  to  a  combination  of  a  linear  mechanism, 
operative  at  low  strains,  and  a  sliding  friction  mechanism, 
producing  nonlinear  effects  at  high  strains.  Thus,  friction 
is  either  masked  by  the  linear  mechanism  or  totally 
inoperative  at  low  amplitudes.  If  the  effect  of  low 
heating  rate  thermal  cycling  is  to  decrease  the 
attenuation  due  to  this  linear  mechanism  and  the  effect 
of  high  heating  rate  cycling  to  increase  attenuation,  then 
the  data  may  be  explained  if  the  contribution  due  to  sliding 
friction  is  constant.  This  mechanism  would  be  appealing 
except  for  the  fact  that  the  linear  mechanism  is  undefined 
in  Mavko' s  model  and  that  cracks,  or  some  grain  contact 
phenomenon  seems  to  be  important  even  at  low  strains.  This 
is  borne  out  by  measurements  made  at  low  amplitudes  (Gardner 
et  al.,  1964;  Winkler  and  Nur,  1978)  where  attenuation 
decreases  with  applied  hydrostatic  pressure  just  as  observed 
in  higher  amplitude  ultrasonic  measurements  (ToksdJz  et  al., 
1979;  Johnston  and  Toksttz,  1979).  All  of  these  data  are 
interpreted  by  the  various  investigators  to  be  due  to  the 
closing  of  thin  cracks,  a  process  that  affects  the  elastic 
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properties  of  rock  as  well  (Walsh,  1965;  Toksfiz  et  al. ,  1976)  . 

An  alternative  explanation  may  be  found  in  variations  of 
crack  contact  area  and  most  easily  thought  of  in  terms  of  work 
done  on  the  contact.  Although  the  small  displacement  at  low 
amplitudes  suggest  that  macroscopic  Coulomb  friction  may  not 
be  strictly  applicable  (Savage,  1969),  it  appears  that  the 
forces  involved  in  attenuation  are  similar.  Walsh  and 
Grosenbaugh  (1978)  have  proposed  a  crack  model  which  holds 
great  promise  to  quantitatively  explain  the  data  reported  in 
this  paper.  The  discussion  here  will  be  qualitative.  Cracks 
are  modeled  as  cavities  formed  by  contact  of  two  bumpy  surfaces. 
The  bumps,  or  asperities,  are  considered  round  and  thus  Hertzian 
contact  theory  is  assumed.  The  distribution  of  asperity 
heights  and  therefore  contact  area  or  normal  stress  is  taken 
to  be  nearly  Gaussian  in  nature.  While  the  exact  form  of  this 
distribution  is  critical  in  terms  of  Walsh  and  Grosenbaugh' s 
elastic  model,  it  is  sufficient  for  this  discussion  that  there 
is  merely  a  peak  in  the  distribution.  A  Gaussian  profile, 
however,  is  observed  for  slightly  polished  steel  and  bead 
blasted  gold  surfaces  (Greenwood  and  Williamson,  1966). 

Considering  the  effect  of  this  model  in  terms  of  general 
amplitude  dependence  of  attenuation,  for  low  strains,  work  is 
done  on  the  low  contact  area  tail  end  of  the  distribution. 

Since  the  distribution  here  is  relatively  flat,  Q  would  appear 


independent  of  amplitude.  As  strain  amplitude  increases,  work 
is  done  on  contact  areas  near  the  peak  of  the  distribution  and 
Q  decreases.  At  very  high  amplitudes,  this  theory  would  predict 
a  leveling  off  of  Q. 

This  same  model  predicts  the  behavior  of  amplitude  effects 
for  thermally  cycled  rocks.  Thermal  cracking  that  widens 
cracks  will  shift  the  peak  in  the  distribution  of  contact  area 
to  lower  values.  This  while  the  actual  density  of  contacts 
may  decrease  and  Q  increases,  the  peak  occurs  at  lower  areas 
and  the  drop'  in  Q  due  to  amplitude  effects  happens  at  lower 
strains . 

The  essential  feature  of  this  model  quite  different  from 
the  Mavko  model,  is  that  contacts  play  a  role  at  all  strain 
amplitude  investigated  in  this  paper.  This  is  more  consistent 
with  our  notion  that  the  increase  in  Q  as  a  function  of  pressure 
measured  at  low  amplitudes  by  resonance  techniques  is  due  to 
contact  closing  just  as  interpreted  for  the  ultrasonic  data. 

All  of  the  amplitude  dependent  models,  however,  require  further 
development. 

In  summary,  the  data  presented  on  tne  variations  of  Q  with 
strain  amplitude  and  thermally  induced  crack  porosity  show  that 
the  attenuation  is  critically  dependent  not  only  on  outgassing 
but  also  the  number  of  cracks  and  their  shapes  or  widths. 

Rocks  thermally  cycled  at  low  heating  and  cooling  rates  up 
to  400 °C  show  a  marked  increase  in  Q  that  cannot  be 


entirely  explained  by  outgassing  of  volatiles.  Such  cycling 
widens  thin  cracks  and  grain  boundaries  thus  reducing  contact 
areas.  In  terms  of  a  sliding  friction-like  mechanism,  this 
would  result  in  lower  attenuation.  Fast  thermal  cycling 
propagates  thin  cracks  and  Q  decreases. 

Q  is  dependent  on  strain  amplitude.  For  low  strains 

_7 

(<10  )  Q  appears  independent  of  amplitude  but  decreases 

—  6  —5 

rapidly  for  strains  greater  than  10  to  10  for  the  rocks 
studied.  There  is  no  hysteresis  in  amplitude  cycling 
suggesting  that  if  friction  is  the  dominant  mechanism  for 
attenuation,  little  or  no  damage  to  the  sliding  surface  occurs 
at  strains  of  10"4  or  smaller.  Amplitude  dependent  behavior 
may  explain  discrepancies  observed  in  Q  values  determined  using 
different  experimental  techniques. 

The  behavior  of  amplitude  dependent  Q  is  also  a  function 
of  the  crack  distribution  in  the  rock.  Slow  thermally  cycled 
samples  experience  the  transition  from  high  to  low  Q  at 
lower  amplitudes  than  virgin  rocks.  Again  this  is  most  easily 
explained  by  appealing  to  variations  in  contact  area  as  a 
result  of  thermal  cycling.  All  of  the  data  imply  that  cracks 
and  grain  boundaries,  and  in  particular,  contact  areas, 
primarily  control  the  attenuation  of  seismic  waves  in  dry  rock. 
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Table  1 


Comparison 

of  Resonance  and 

Ultrasonic 

Q  Data 

Berea 

Navajo 

Resonance, 

low  amplitude 

75 

60 

Resonance , 

high  amplitude 

30 

28 

Ultrasonic 

P  wave 

20 

30 

Table  2 

Crack  Porosity  (%)  From  Thermal  Cycling 


Virgin  400°C  300°C 
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Solenhofen  Limestone  —  0.09  0.18 

Westerly  Granite 

perpendicular  0.30  0.54  2.70 

parallel  0.12  0.34  2.07 
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FIGURE  CAPTIONS 


1.  Q  and  bar  velocity  as  a  function  of  strain 
amplitude  for  longitudinal  waves  in  the  dry  Berea 
sandstone.  Also  shown  is  Q  as  a  function  of  strain 
amplitude  for  plexiglass.  In  this  and  later  figures, 
representative  error  bars  are  shown. 

2.  Q  as  a  function  of  strain  amplitude  for  longitudinal 
waves  in  the  dry  Navajo  sandstone. 

3.  Linear  strain  as  a  function  of  confining  pressure  in 
Frederick  diabase  for  a  virgin  sample  and  samples 
thermally  cycled  to  400  and  800°C. 
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Solenhofen  limestone  for  virgin  and  thermally  cycled 
samples. 

Fig.  5.  Linear  strain  as  a  function  of  confining  pressure  in 
Westerly  granite  for  the  "soft"  direction.  Curves  for 
the  virgin  and  400°C  cycled  samples  are  shown. 

Fig.  6.  Linear  strain  as  a  function  of  confining  pressure  in 
Westerly  granite  for  the  "stiff"  direction. 

Fig.  7.  Linear  strain  as  a  function  of  confining  pressure  in 
Westerly  granite  for  both  directions  for  samples 
thermally  cycled  to  800°C. 


Fig.  8.  Q  (solid  symbols)  and  bar  velocity  (open  symbols) 
as  a  function  of  maximum  temperature  achieved  in  slow 
thermal  cycling  in  the  Frederick  diabase.  Two  separate 
experiments  were  run  and  the  results  are  denoted  by  the 
different  symbols.  Representative  error  bars  are  also 
shown . 

Fig.  9.  Q  and  bar  velocity  as  a  function  of  maximum 

temperature  in  Solenhofen  limestone.  See  Figure  8  for 
explanations . 

Fig.  10.  Q  and  bar  velocity  as  a  function  of  maximum 

temperature  in  Westerly  granite  for  both  directions. 

See  Figure  8  for  explanations. 

rig.  il.  Q  (soxxd  symooxs;  ana  oar  vexocity  (open  symooxs; 
as  a  function  of  maximum  temperature  achieved  in  fast 
thermal  cycling  in  Solenhofen  limestone. 

Fig.  12.  Q  and  bar  velocity  as  a  function  of  maximum 
temperature  (fast  cycling)  in  Frederick  diabase. 

Fig.  13.  Q  as  a  function  of  strain  amplitude  in  Frederick 
diabase  for  virgin  and  several  slow  thermally  cycled 
samples . 

Fig.  14.  Q  as  a  function  of  strain  amplitude  in  Solenhofen 
limestone  for  virgin  and  several  slow  thermally  cycled 
samples. 


Fig.  15.  Q  as  a  function  of  strain  amplitude  in  Westerly 
granite  (soft  direction)  for  virgin  and  several  slow 


thermally  cycled  samples. 

16.  Comparison  of  amplitude  dependent  Q  behavior  for 
slow  and  fast  cycled  samples  of  Frederick  diabase. 
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ABSTRACT 


Ultrasonic  P  and  S  wave  Q  values  as  functions  of 
hydrostatic  pressure  are  presented  for  the  Berea  and  Navajo 
sandstones,  Bedford  limestone,  two  shales  and  two  tuffs. 

In  most  cases,  the  attenuation  was  obtained  for  both  dry 
and  water  saturated  samples.  Q  is  determined  by  measuring 
the  attenuation  of  the  rock  sample  relative  to  a  high  Q 
standard  of  identical  geometry  using  Fourier  spectral 
ratios.  The  data  show  that  Q  increases  with  pressure, 
the  rate  of  increase  being  dependent  on  rock  type  and 
crack  density  and  distribution.  Qs  is  about  equal  to  or 
larger  than  C,~  for  dry  rocks  and  less  than  Q,>  for  saturated 
rocks.  In  all  cases,  Q's  for  saturated  rocks  are  lower 
than  for  dry  samples.  The  variations  of  attenuation  in 
shales  are  not  crack  controlled  and  an  oil  shale  studied 
shows  extreme  and  unusual  anisotropy  in  velocity  and  Q. 

At  ultrasonic  frequencies ,  the  attenuation  in  the  tuffs 
studied  is  dominated  by  scattering  losses. 


INTRODUCTION 


The  attenuation  of  seismic  waves  has  great  potential 
as  a  tool  to  understand  the  physical  conditions  of  rocks  in 
the  upper  crust.  Unfortunately,  the  systematics  of  the 
behavior  of  attenuation  with  pressure  and  saturation 
conditions  are  not  well  studied.  In  this  paper  we  present 
'  ultrasonic  Q  data  as  functions  of  pressure  for  several  dry 
and  saturated  crustal  rocks. 

Attenuation  measurements  have  been  carried  out  on  samples 
in  the  laboratory  using  a  variety  of  techniques  (Attewell  and 
Ramana,  1966;  Birch  and  Bancroft,  1938;  Born,  1941;  Bradley 
and  Fort.  19C»6;  Car-Hn^r-  et  al.,  1964;  Hamilton.  1972;  Jackson. 
1969;  Knopoff,  1964;  Nur  and  Simmons,  1969;  Peselnick  and 
Zietz,  1959;  Peselnick  and  Outerbridge,  1961;  Spetzler  and 
Anderson,  1968;  Tittmann  et  al.,  1974;  Toksfiz,  et  al.,  1979; 
Warren  et  al.,  1974;  Watson  and  Wuenschel,  1973;  Wyllie  et 
al.,  1962).  In  a  most  general  way,  these  measurements  show 
that  Q  is  independent  of  frequency,  fluid  saturation  increases 
attenuation  (Gardner  et  al.,  1964;  Obert  et  al.,  1946;  Watson 
and  Wuenschel,  1973;  Wyllie  et  al.,  1962),  and  increasing 
pressure  decreases  attenuation  (Gardner  et  al.,  1964;  Gordon 
and  Davis,  1968;  Klima  et  al.,  1964;  Levykin,  1965;  ToksOz 
et  al.,  1979).  However,  more  laboratory  data  are  needed 
under  controlled  conditions  to  determine  the  effects  of 
fluids  and  pressure  on  the  attenuation  in  porous  rocks. 


EXPERIMENTAL  METHOD 


Q  values  for  P  and  S  body  waves  are  determined  with  an 
ultrasonic  pulse  transmission  technique  described  in  Toksfiz 
et  al.  (1979) .  The  attenuation  is  measured  relative  to  a 
reference  sample  of  low  attenuation  using  spectral  ratios. 
Pulse  transmission  methods  are  superior  to  resonant  bars  for 
measurements  taken  on  jacketed  samples  under  pressure  if 
geometric  factors  such  as  beam  spreading  and  reflections 
can  be  corrected  for.  Only  one-way  transmission  effects 
are  measured.  The  sample  to  be  studied  arid  reference  sample 
have  exactly  the  same  shape  and  geometry.  Essentially,  two 
measurements  are  made  using  identical  procedures,  one  with 
me  rocx  or  mceresC  -  ciiiu  cigaxn  wxui  Liic 

The  amplitudes  of  plane  seismic  waves  for  the  reference 


and  the  sample  can  be  expressed  as 


i  /  r  \  r  ,  \  ) 

A,  ($)  =  6,(x)e  e 
Aj.(f')  --  G,  U) e'^'A 


where  A  =  amplitude,  f  =  frequency,  x  =  distance,  k  =  2Uf/v 
=  wavenumber,  v  =  velocity,  and  G(x)  is  a  geometrical  factor 
which  includes  spreading,  reflections,  etc.  <*(f)  is  the 
frequency  dependent  attenuation  coefficient.  Subscripts  1 
and  2  refer  to  the  reference  and  sample  respectively.  From 


available  data  it  is  reasonable  to  assume  that  over  the 
frequency  range  of  measurements,  0. 1-1.0  MHz,  «  is  a  linear 
function  of  frequency,  although  the  method  itself  tests  this 
assumption  {Knopoff,  1964;  Jackson  and  Anderson,  1970; 
McDonal,  et  al.,  1958).  Then  one  can  write; 

*(0  -  y-f 


where  is  constant  and  related  to  the  quality  factor  Q  by 


When  the  same  geometry  is  used  for  both  the  sample  and  standard 
(i.e. ,  same  sample  dimensions,  transducers,  arrangements)  then 
and  G2  are  frequency  independent  scale  factors.  The  ratio 
of  the  Fourier  amplitudes  are; 


A . 
Az 

A. 

A* 


6,  (^1  ^  x 
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(y»  -  Ir.Hx  +  ^  (f) 


Since  G^A»2  is  independent  of  frequency  and  x  is  the  sample 


f 


129 


length,  (Y^-  83.)  can  be  foun<^  from  the  slope  of  the  line  fitted 
to  lnCA^/^)  versus  frequency.  If  the  Q  of  the  standard 
reference  is  known,  ^2  of  the  sample  can  be  determined.  When 
the  Q  of  the  standard  is  very  high  (i.e.  =«*►)  then 

Y^  =  0  and  >^2  the  rock  sample  can  be  determined  directly 
from  the  slope. 

The  samples  used  in  this  study  are  cylindrical,  2.54  cm 
long  and  3.175  cm  in  diameter.  Transmitter  and  receiver 
transducers  (each  2.5  cm  in  diameter)  are  mounted  in  a  stacked 
P  and  S  wave  arrangement  at  opposite  ends  of  jacketed  samples. 
The  pore  or  fluid  pressure,  P^,  system,  capable  of  2  kb, 
is  controlled  independently  of  the  confining  pressure,  Pc 
(rated  to  10  kb).  Generally,  the  experiments  are  run  to  a 
maximum  differential  pressure,  Pd  =  Pc  -  P^,  of  2.25  kb  in 
discrete  steps.  At  each  pressure  step,  velocity  for  the 
ultrasonic  pulse  is  determined  by  matching  an  electronically 
delayed  timing  pulse  with  the  first  arrival,  transit  time 
being  measured  to  a  10  ns  accuracy  using  a  time  interval 
counter.  Full  P  and  S  waveforms  are  recorded  at  each  step 
and  then  manually  digitized.  These  waveforms  are  then 
windowed  and  cosine  tapered  at  the  end  to  remove  reflections. 

A  standard  FFT  routine  is  used  to  obtain  the  amplitude  spectra. 

In  this  study  we  used  aluminum  as  a  standard  reference. 

Q  for  aluminum  is  about  150,000  (Zamanek  and  Rudnick,  1961), 
as  opposed  to  Q  £  1000  for  rocks.  Thus  taking  =  0  never 


130 


r~ 


-r ,  I-; 


v.v;, 


introduces  more  than  1%  error.  For  typical  rocks  where 
Q  =  10-100,  the  error  is  less  than  0.1%  and  is  negligible. 

A  more  serious  concern  is  the  validity  of  assumption  that  the 
geometric  factors  and  G2  have  the  same  frequency 
dependence,  and  G^/G2  is  independent  of  frequency.  With 
polished  rock  surfaces  and  good  coupling  between  the 
transducer  holder  and  sample,  one  would  not  expect  frequency 
dependent  reflection  coefficients  at  the  interface. 

Formal  errors  of  the  Q  measurements,  however,  are 
difficult  to  assign.  Strict  interpretation  of  the  errors  on 
fitting  a  straight  line  to  the  amplitude  ratios  only  determines 
the  modulation  character  of  the  ratios  and  not  the  accuracy  of 
the  method.  It  is  clear,  however,  that  a  reliable  estimate  or 
Q  may  be  obtained  only  if  the  slope  of  ln(A1/A2)  versus 
frequency  is  large.  Thus  the  method  is  not  satisfactory  for 
samples  with  high  Q  (>100) . 

The  resolution  and  accuracy  of  the  method  can  be  based 
only  on  reproducibility  and  comparison  with  other  techniques. 
Repeated  measurements  showed  that  pulse  amplitudes,  shapes 
and  spectra  were  duplicates  and  Q  values  determined  to  within 
about  15%.  Digitizing  errors  were  found  to  be  no  larger  than 
the  repeatability.  Furthermore,  results  described  in  this 
paper  for  the  Berea  sandstone  may  be  compared  to  data  (Toksdz 
et  al.,  1979)  for  the  same  rock  at  the  same  pressure  and 
saturation  conditions  using  the  same  technique  but  obtained  at 
a  different  laboratory  with  different  sample  dimensions. 


•v.v.n'.v.v 
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There  is  excellent  agreement.  Thus  the  reproducibility  of  the 

f 

Q  measurements  may  be  conservatively  estimated  at  15-20%.  To 
check  for  bias  in  the  method,  Q  was  obtained  for  plexiglas 
using  the  ultrasonic  pulse  and  resonant  bar  techiques.  The 
Q  values  determined,  20  and  22  respectively,  agreed  to  within 
the  resolution  of  the  spectral  ratio  method  and  thus  we  may 
assume  are  not  overly  biased.  Error  bars  for  several 
reprentative  data  points  are  shown  on  each  graph  in  this 
paper. 

As  examples  of  the  Q  measurement  technique ,  waveforms , 
normalized  spectra  and  spectral  ratios  are  shown  for  P  and  S 
waves  in  the  dry  sandstone  in  Figures  la.b  for  two 

confining  pressures.  In  these  figures,  the' dashed  line  in 
the  spectrum  plots  corresponds  to  the  aluminum  standard  (Al) . 
Only  waveforms  for  the  rock  sample  are  shown.  A  rather  large 
DC  shift  is  noted  for  the  S  waves.  This  feature,  common  to 
many  of  the  experiments,  is  due  to  low  frequency  noise.  This 
could  be  filtered  out  prior  to  time  windowing  but  was  not 
necessary  since  the  low  frequency  portion  of  the  spectrum  is 
not  used  in  determining  Q.  One  can  see  from  the  figures, 
though,  that  the  slopes  of  the  spectral  ratio  plots  clearly 
decrease  with  increasing  pressure  indicating  an  increase  in 
Q.  Q  values  are  calculated  from  these  slopes  and  the  measured 
P  and  S  wave  velocities. 


DESCRIPTION  OF  SAMPLES 

The  rock  samples  chosen  for  this  study  include  two 
sandstones,  a  limestone,  two  shales  and  two  tuffs.  Relatively 
low  Q's  were  expected  in  these  rocks  making  them  amenable  to 
the  pulse  technique.  Furthermore,  most  of  the  samples  are 
quite  commonly  used  in  the  laboratory  and  thus  much  subsidiary 
information  on  their  properties  may  be  obtained.  In  fact, 
the  Berea  sandstone  and  Bedford  limestone  are  accepted 
"standard"  rocks  by  the  U.S.  Bureau  of  Mines  (Krech  et  al., 
1974). 

Each  sample  was  characterized  in  terms  of  mineral  content, 
grsin  size  and  contacts,  density,  porcsitv,  permeability  and 
static  stress-strain  properties.  Detailed  sample  descriptions 
may  be  found  in  Johnston  (1978) .  Several  important  physical 
parameters  for  each  rock  are  listed  in  Table  1.  Specifically 
the  following  rocks  were  studied: 

1.  Berea  Sandstone.  The  Berea  is  a  relatively  clean, 
medium  to  fine-grained  protoquartzite  with  silica  and  clay 
cement.  The  porosity  of  the  sample  used  was  18%  with  a 
permeability  of  about  75  md.  Its  static  response,  typical  of 
many  rocks,  is  shown  in  Figure  2  where  linear  strain  is  plotted 
as  a  function  of  hydrostatic  confining  pressure.  The  non-linear 
portion  of  the  curve  is  due  to  crack  closure  (Walsh,  1965)  and 
permanent  hysteresis  upon  unloading  indicates  that  not  all  the 
cracks  have  closed  elastically. 

2.  Navajo  Sandstone.  The  Navajo  is  classified  as  an 
orthoquartzite  cemented  primarily  by  silica  with  only  minor 


amounts  of  clay.  The  sample  studied  had  a  porosity  of  16% 
and  an  estimated  permeability  of  100  md.  The  static  response 
is  very  similar  to  that  observed  for  the  Berea  sandstone. 

3.  Bedford  Limestone.  Also  known  as  the  Salem,  Spergen 
and  Indiana  limestone,  the  Bedford  is  coarse  grained,  made  of 
poorly  sorted  fossil  debris  cemented  by  crystalline  calcite. 

The  porosity  of  the  sample  used  was  12%  and  the  permeability 
is  about  1  md  (Krech  et  al.,  1974).  Pore  collapse  under 
hydrostatic  pressure  occurs  at  about  1  kb  as  shown  in  Figure  3. 
Upon  unloading  the  sample,  a  large  increase  in  crack  porosity 
due  to  the  pore  collapse  is  noted.  It  will  be  seen  that  this 
has  a  large  effect  on  the  acoustic  properties. 

4.  Esopus  Shale.  The  Esopus  is  a  dark  gray,  calcareous 
silty  shale  found  in  eastern  New  York  state.  The  porosity  and 
permeability  were  both  too  small  to  be  measured.  Furthermore, 
when  dry,  the  rock  is  extremely  friable  and  thus  only  a 
saturated  experiment  was  run. 

5.  Colorado  Oil  Shale.  The  sample  used  in  this  study 

was  obtained  from  the  Anvil  Points  Mine  and  represents  a  typical 
medium  grade  oil  shale  with  about  20%  kerogen  by  weight.  Again 
there  is  no  measureable  porosity  or  permeability.  The  rock 
is  visually  aniostropic  and  cores  were  taken  both  normal  and 
perpendicular  to  the  bedding  planes.  Static  measurements  show 
a  30%  anisotropy  in  linear  strain.  Little  or  no  crack  porosity 
is  evident  and  there  is  little  permanent  hysteresis. 


6.  Ammonia  Tanks  Tuff.  The  Ammonia  Tanks  is  classified 
as  a  moderately  welded  quartz  latite  and  is  compositionally 
complex  (Byers  et  al.,  1976).  The  porosity  of  the  sample  used 
was  measured  as  5.8%  but  this  probably  does  not  include 
substantial  unconnected  porosity.  Grain  sizes  can  be  as  large 
as  5  mm.  The  static  response  of  this  and  the  following  tuff 
at  low  stresses  (4  60  bars)  suggest  that  compared  to  the 
sandstones,  far  fewer  fine  cracks  exist  (Morrow  et  al.,  1977). 
The  total  crack  porosity,  however,  is  larger  (Johnston,  1978) 
implying  that  the  bulk  of  the  cracks  in  the  tuffs  are  fairly 
wide. 

7.  Tunnel  Beds  Tuff.  The  Tunnel  Beds  is  an  ash  fall  tuff, 
poorly  sorted  with  very  coarse  grains  2-4  mm  in  diameter.  The 
measured  porosity  was  29.6%.  As  with  the  Bedford  limestone, 
pore  collapse  occurs  under  hydrostatic  compression  but  at  a 
pressure  of  about  500  bars. 

EXPERIMENTAL  DATA 

In  this  section,  ultrasonic  Q  data  as  functions  of 
saturation  and  pressure  are  presented.  For  each  experiment 
an  aluminum  sample  with  the  same  shape  of  the  rock  sample 
was  used  as  a  standard  reference.  Q  values  were  calculated 
from  spectral  ratios  of  rock  -  aluminum  pairs  at  each  pressure. 

1 .  Berea  Sandstone 

The  Berea  sandstone  is  a  particularly  well-studied  rock 
and  Q  values  as  functions  of  pressure  and  saturation  have  been 
published  in  Toksfiz  et  al.  (1979).  Here,  those 


experiments  are  repeated  and  several  new  ones  run.  In  all, 
four  experiments  were  run  on  the  Berea:  one  with  a  dry  sample, 
the  second  fully  water  saturated,  the  third  partially  saturated 


increases.  However,  the  trend  in  Q  versus  pressure  is  clear. 

Qc  is  about  equal  to  Q  . 

P 

A  Berea  sample  was  then  saturated  with  distilled  water 
and  the  experiment  repeated.  The  fluid  pressure  was  maintained 
at  0.465  the  confining  pressure  for  differential  pressures 
ranging  from  45  to  2000  bars.  This  ratio  of  pore-fluid  and 
confining  pressures  is  a  nominal  value  for  saturated  sedimentary 
rocks  where  the  fluid  pressure  is  equal  to  the  hydrostatic 
pressure  of  a  water  column. 

The  presence  of  water  increases  the  velocity  of  P  waves 

while  decreasing  the  velocity  of  S  waves  as  illustrated  in 

Plgnra  4.  The  increase  in  velccitv  with  pressure  is  not  as 

great  as  for  the  dry  sample.  Q  values  are  shown  in  Figure  6 

as  functions  of  differential  pressure.  Both  Qp  and  Qg  are 

substantially  lower  than  their  dry  counterparts  and  in  the 

saturated  case  Q_  is  much  lower  than  Q_,  approaching  a  value 

5  P 

of  only  about  30  even  at  high  pressure.  Qp  increases  to  a 
value  of  about  60  by  a  pressure  of  about  750  to  1000  bars. 

While  it  appears  that  the  effect  of  fluid  saturation  is  to 
dramatically  increase  the  attenuation,  especially  for  S  waves, 
the  dependency  of  the  Q  values  on  pressure  is  still  observed. 

A  similar  set  of  Q  determinations  for  the  Berea  sandstone 
were  obtained  by  Toksflz  et  al.  (1979).  While  the  physical 
conditions  (pressure  and  saturation)  of  the  experiments  were 
similar  to  those  presented  here,  the  geometry  was  somewhat 
different.  The  results  of  the  two  sets  of  experiments  agree 
very  well,  enhancing  our  confidence  in  the  experimental  technique. 


Finally/  two  experiments  were  run  to  investigate  the 


effects  of  partial  water  saturation  on  Q.  The  primary 
influence  of  water  saturation  in  known  from  data  like  those 
presented  by  Gardner  et  al.  (1964)  for  the  Berea.  Q  for  both 
longitudinal  and  torsional  vibrations  of  bars  was  observed  to 
decrease  rapidly  relative  to  a  dry  sample  for  small  amounts  of 


water  saturation.  In  these  experiments  we  start  from  fully 
saturated  samples  and  then  reduce  the  water  saturation. 

In  the  first  run,  partial  saturation  was  obtained  by 
partially  drying  out  a  completely  saturated  sample.  The  sample 


was  stored  in  an  airtight  plastic  bag  for  over  a  week  to  allow 


jacketing,  the  sample  was  reweighed  and  the  saturation  was 
determined  to  be  37%.  Q  values  for  this  experiment  are  shown 
in  Figure  7.  Comparing  these  data  with  those  of  the  completely 
water  saturated  Berea  sandstone  (Figure  6) ,  one  finds  that  Qp 
is  slightly  lower  relative  to  the  fully  saturated  sample  at 
intermediate  pressures  while  Qs  is  essentially  unchanged. 

There  is  no  evidence,  however,  that  P  wave  attenuation  changes 
more  than  a  factor  of  two  relative  to  full  saturation  at  this 
degree  of  partial  saturation. 

In  order  to  examine  this  in  greater  detail,  an  experiment 
was  designed  so  that  partial  saturation  could  be  achieved  by 
the  release  of  helium  gas  from  saturated  pore  water  by 
reducing  the  fluid  pressure  in  a  closed  system.  The  pore 
pressure  apparatus  and  the  procedure  used  for  this  experiment 


is  described  in  Johnston  (1978).  At  a  fluid  pressure  of 
200  bars  and  a  confining  pressure  of  400  bars,  the  completely 
water  saturated  sample  and  pore  pressure  system  was  injected 
with  the  volume  of  gas  required  to  completely  saturate  the 
water  at  a  fluid  pressure  of  1000  bars.  Fluid  pressure  was 
then  increased  to  1200  bars  maintaining  a  differential 
pressure  of  200  bars  by  simultaneously  raising  the  confining 
pressure  to  1400  bars.  Then,  fluid  and  confining  pressure 
were  both  lowered,  allowing  the  gas  to  come  out  of  solution. 
Starting  from  zero  gas  saturation  at  1000  bars  fluid  pressure, 
a  total  gas  saturation  of  5%  was  obtained  at  a  fluid  pressure 
of  150  bars.  Velocities  were  determined  and  waveforms  recorded 
at  discrete  pressures. 

The  presence  of  free  gas  in  the  rock  was  confirmed  two 
ways.  First,  the  compressibility  of  the  pore  fluid  could  be 
roughly  determined  by  the  volume  change  required  to  affect  a 
given  pressure  drop.  This  clearly  indi  rted  an  increase  in 
compressibility  at  fluid  pressures  below  1000  bars  but  assumed 
that  the  gas  phase  was  homogeneously  distributed  throughout 
the  system.  More  concrete  evidence  was  found  in  the  velocities 
shown  in  Figure  8.  Here,  and  Vg  are  plotted  as  functions 
of  pore  fluid  pressure  and  calculated  gas  saturation.  Clearly, 
Vp  decreases  with  increasing  gas  saturation  while  Vg  remains 
constant.  This  behavior  is  exactly  as  would  be  expected  from 
theory  but  as  a  check,  Vp  was  determined  theoretically 


(Kuster  and  Toks6z,  1974;  Toksfiz  et  al. ,  1976)  for  a  change 
was  the  compressibility  of  the  water  as  a  function  of  fluid 
pressure.  This  calculation  is  shown  as  the  dashed  line  in 
Figure  8  and  obviously  does  not  explain  the  data  except  for 
the  two  high  pressure  points  which  are  supposed  to  be  completely 
water  saturated  anyway.  It  appears  then,  that  the  experiment 
was  indeed  successful  in  releasing  helium  gas  from  solution 
in  the  rock. 

Q_  and  Q  are  shown  as  functions  of  gas  saturation  and 
p  s» 

fluid  pressure  in  Figure  9.  The  average  values  of  both  are 
comparable  to  those  obtained  at  200  bars  differential  pressure 
in  the  previous  experiments.  Clearly,  there  is  no  change  in 
Qg  as  a  function  of  gas  saturation  in  0  to  5%  range.  There  are 

t 

variations  in  Q  and  even  though  they  are  within  the  resolution 

r 

of  the  experimental  technique,  they  might  be  significant.  In 
particular,  there  is  a  reduction  in  Qp  up  to  a  gas  saturation 
of  3  or  4%  at  which  Qp  increases  to  its  full  water  saturation 
value . 

While  the  significance  of  these  results  are  questionable, 
it  should  be  pointed  out  that  the  maximum  effects  of  partial 
gas  saturation  are  to  be  expected  at  low  pore  pressures.  Thus, 
the  changes  in  Q  observed  here  might  possibly  be  amplified 


somewhat. 


2.  Navajo  Sandstone 


As  with  the  Berea  sandstone,  Q  values  were  obtained  for 
the  Navajo  in  a  vacuum  dried  sample  and  a  water  saturated 
sample.  Velocity  as  a  function  of  pressure  is  shown  in 
Figure  10  for  P  and  S  waves.  Consistent  with  static  stress- 
strain  measurements,  velocity  increases  more  rapidly  at  lower 
pressures  than  for  the  Berea.  Furthermore,  the  values  for 
velocity  are  higher  -  3.56  km/s  versus  3.22  km/s  for  P  waves 
in  dry  rock  at  45  bars  and  4.37  km/s  versus  4.23  km/s  at 
1000  bars  as  examples.  Similar  comparisons  may  be  made  for 
the  saturated  samples  and  for  S  waves. 

Q  vclucs  f cr  t.hc  drv  rcch  cbtsin^d  f— on1,  cnorf rai 

ratios  are  shown  in  Figure  11.  Again,  it  must  be  mentioned 
that  the  high  Q  values  at  the  higher  pressures  are  unreliable 
a  restriction  of  the  experimental  technique.  However,  for 
pressures  less  than  500  bars,  one  may  say  that  the  Q’s 
calculated  for  both  P  and  S  waves  behave  in  much  the  same 
manner  as  for  the  dry  Berea  sandstone.  The  absolute  values 
of  Q  are  higher  than  for  the  Berea  -  30  compared  to  24  for  P 
waves  at  45  bars  confining  pressure.  Q  also  seems  to  increase 
at  a  slightly  faster  rate  with  pressure. 

In  the  second  experiment  the  sample  was  saturated  with 
distilled  water  and  as  with  one  of  the  Berea  runs,  P^  was 
maintained  at  0.465  Pc.  Q  values  are  shown  in  Figure  12. 

In  this  case,  the  differences  between  the  Navajo  and  Berea 
sandstones  are  more  clearly  observed.  While  the  Q's  obtained 


for  the  Navajo  are  substantially  the  same  as  those  for  the 
Berea  at  pressures  greater  than  600-700  bars,  the  increase  in 
Q  at  low.  pressures  is  much  faster  for  the  Navajo.  Also,  as 
with  the  dry  case,  the  absolute  Q  values,  particularly  for 
P  waves,  are  higher  for  the  Navajo.  This  is  in  good  agreement 
with  the  differences  observed  in  the  velocity  measurements. 
Clearly,  this  is  a  good  example  of  how  the  same  parameters 
that  determine  elastic  properties  also  affect  attenuation. 

3.  Bedford  Limestone 

The  Bedford  limestone  presents  an  interesting  comparison 
to  the  sandstone.  Although  the  total  porosity  is  rather  high 
(12%),  the  crack  porosity  is  somewhat  lower  than  the  sandstone. 
AioO ,  pore  collapcc,  initiated  at  I'J'J'J  tare.  ? • 
behavior  of  both  velocity  and  Q  observed  in  the  sandstones. 
Again,  two  experiments  were  run,  one  dry  and  the  other 
saturated.  For  the  water  saturated  case  Pf  was  maintained  at 
.465  Pc  throughout  the  experiment. 

P  and  S  wave  velocities  for  the  Bedford  limestone  are 
shown  in  Figure  13.  We  note  the  general  features  of  an 
increase  in  velocity  due  to  crack  closure  up  to  about  500  bars 
and  higher  P  wave  and  lower  S  wave  velocities  for  the  saturated 
sample.  At  1000  bars  the  effect  of  pore  collpase  becomes 
evident.  During  the  loading  cycle  (increasing  differential 
pressure) ,  velocity  increases  as  the  large  pores  collapse. 

As  the  rock  is  unloaded  (decreasing  pressure) ,  it  behaves 
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elastically  until  about  800  bars  at  which  point,  cracks 
induced  by  pore  collapse  open,  reducing  the  velocities  to 
a  point  well  below  those  obtained  during  the  loading  cycle. 
This  type  of  behavior  was  identified  by  Nur  (1969)  and,  of 
course,  confirmation  lies  in  the  static  stress-strain 
measurements  shown  in  Figure  3.  The  effect  is  observed  for 
both  P  and  S  waves  in  both  dry  and  saturated  samples. 

There  was  some  question  as  to  the  quality  of  the  shear 
waves  in  the  dry  rock  and  thus  the  results  should  be  used 
with  caution.  As  we  will  see,  the  dry  Qs  values  are  not 
consistent  with  the  other  data.  There  is  no  explanation  at 
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values  for  the  dry  sample  are  plotted  as  functions  of  confining 
pressure  (Pf  =  0)  in  Figure  14.  Concentrating  on  the  P  wave 
values,  we  see  that  Qp  is  much  higher  than  for  sandstones  at 
low  pressures  -  around  70  at  105  bars.  The  rate  of  increase 
with  pressure  is  lower  than  the  sandstones  however,  and  at 
pressures  of  700-800  bars,  Q  is  only  slightly  higher  than  the 

r 

values  attained  by  the  sandstones.  While  the  data  is  shown 
only  to  a  pressure  of  1000  bars  (subjected  to  the  constraints 
of  the  experimental  method) ,  the  sample  was  run  to  a  maximum 
pressure  of  2250  bars.  During  the  unloading  cycle,  the  effect 
of  pore  collapse  is  evident  at  the  lower  pressures.  Qp  is 
sharply  reduced  at  pressures  less  than  1000  bars  due  to  the 
large  increase  in  crack  porosity.  For  example,  at  100  bars 
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the  Q  is  lowered  by  52%  (118%  increase  in  attenuation) .  This 
may  be  compared  to  the  344%  increase  in  crack  porosity 
obtained  from  the  static  measurements.  Clearly,  not  all  of 
these  cracks  are  of  sufficiently  low  aspect  ratio  (crack 
width  to  length  ratio)  at  100  bars  to  contribute  to  the 
attenuation.  Indeed,  Q  appears  to  level  off  from  250  to  100 
bars  on  the  unloading  cycle  implying  a  lower  limit  to  the 
induced  crack  aspect  ratio  spectrum  at  surface  pressure. 

Q  values  calculated  for  the  saturated  Bedford  limestone 

are  shown  in  Figure  15.  Again  the  Q  is  lower  than  for  the 

dry  case  -  by  a  factor  of  two  for  P  waves  at  500  bars  for 

example,  and  Q  is  lower  than  Q  .  However,  we  see  little 
-  P 

difference  between  the  supposed  dry  Qs  values  (Figure  14) 
and  the  wet  ones.  Nor  do  the  dry  S  results  increase  with 
pressure  as  the  wet  data.  We  would  expect  from  what  we  have 
seen  in  the  velocities,  that  the  changes  in  Q  for  the  dry 
rock  at  pressures  less  than  500  bars  would  be  more  rapid 
than  the  wet  rock. 

In  any  event,  the  rapid  increase  in  Q  at  low  pressure 
is  followed  by  a  levelling  off  up  to  about  900  bars.  An 
increase  in  t,  particularly  for  S  waves,  occurs  at  that  point 
due  to  pore  collapse.  The  final  value  of  Qg  obtained  at  2000 
bars  is  nearly  twice  that  obtained  by  the  sandstones.  The 
increase  in  Qp  is  not  as  dramatic.  However,  upon  unloading 
the  sample,  the  effect  of  induced  cracks  is  much  greater  for 
P  waves  than  S  waves  although  both  obtain  lower  values  than 
during  the  loading  cycle. 


4.  Shales 


The  first  experiment  run  was  on  the  Esopus  Shale.  Because 
the  porosity  and  permeability  of  this  sample  is  so  low,  a  full 
series  of  pressures  was  not  used.  Furthermore,  since  the 
rock  becomes  very  friable  when  dry,  only  a  water  saturated 
test  was  run.  Rather  than  applying  pore  pressure,  the  sample 
remained  with  atmospheric  pore  pressure  during  the  course  of 
the  experiment. 

Data  were  obtained  at  105,  500,  1000,  and  2000  bars 
confining  pressure.  The  sample  was  held  at  each  pressure 
for  two  hours  to  allow  the  pore  pressure  to  equilibrate. 

Vr  1  r.r”i  tins  vr.ro  continuously  monito^®^  dnrH  nn  this  time 
and  the  experiment  continued  once  the  velocity  leveled  off 
to  a  constant  value.  However,  the  sample  was  kept  at  1000  bars 
overnight  (16  hours)  resulting  in  a  further  increase  in  velocity. 
This  implied  that  2  hours  were  not  enough  to  equilibrate 
pore  pressure.  Velocities  obtained  unloading  the  sample  show 
that  the  data  from  1000  bars  (after  the  16  hours  wait)  and 
2000  bars  probably  most  closely  represent  the  completely 
equilibrated  values. 

V_  and  V_  for  the  2  hour  equilibrated  points  and  the 

hr  » 

16  hour  equilibrated  run  at  1000  bars  are  shown  as  functions 
of  confining  pressure  in  Figure  16.  The  changes  in  velocity 
are  not  nearly  as  .great  as  for  the  sandstones  or  limestones 
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due  to  the  absence  of  fine  cracks  implied  by  the  lack  of 
effective  porosity  or  permeability.  The  unloading  portions 
of  the  velocity  curves,  while  offset  from  the  low  pressure 
loading  curves  due  to  pore  pressure  equilibration,  have 
essentially  the  same  curvature,  further  corroborating  the 
absence  of  fine  cracks. 

Very  little  obvious  change  in  the  nature  of  the  waveforms 
and  spectra  was  observed  implying  that  Q  does  not  change  very 
drastically.  This  is,  in  fact,  clear  in  Figure  17  where  the 
Q  values  obtained  from  the  spectral  ratios  are  plotted  as 
functions  of  confining  pressure.  The  effect  of  pore  pressure 

t 
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velocity.  The  Q  values  for  the  Esopus  shale  are  consistently 

higher  than  for  the  saturated  sandstones  at  all  pressures 

although  the  difference  is  minimized  at  higher  pressures  where  * 

the  fine  cracks  in  the  sandstones  have  closed.  The  large 

difference  between  Q  and  Q  may  be  due  to  slippage  in  clay 

P  s 

lubricated  grain  contacts  oriented  parallel  to  the  bedding 
and  perpendicular  to  the  direction  of  S  wave  propagation. 

Both  Qp  and  Qs  are  also  higher  than  in  situ  values  (32  and  10 
respectively)  obtained  for  the  Pierre  shale  by  McDonal  et  al. 
(1958)  but  lie  well  within  the  range  of  Q  values  for  various 
shales  listed  by  Bradley  and  Fort  (1966) . 

Two  experiments  were  run  for  the  Colorado  Oil  Shale  - 
one  perpendicular  to  the  bedding  and  the  other  parallel. 


As  with  the  Esopus  shale  the  samples  were  completely  saturated 
but  P^  =  0  throughout  the  experiment.  The  velocities,  plotted 
in  Figure  18  are  consistent  with  the  static  measurements 
mentioned  earlier.  For  both  P  and  S  waves,  velocity  is 
higher  in  the  direction  parallel  to  the  bedding.  There  is 
little  evidence  of  microcracks  since  the  velocities  show  no 
rapid  changes  at  low  pressures  and  increase  less  than  20% 
by  2000  bars  compared  to  an  increase  of  over  30%  for  sandstones. 

Because  the  slopes  of  velocity  versus  pressure  are  virtually 
identical  for  both  orientations ,  we  may  assume  that  the 
variations  in  velocity  observed  are  not  controlled  by  cracks. 
Furthermore,  the  presence  of  krrogen  in  the  shale  seems  to 
play  an  important  role  in  determining  both  the  elastic  and 
anelastic  properties.  First,  the  Poisson's  ratio  calculated 
from  the  velocities  is  0.40  -  implying  that  the  shale  is 
behaving  more  like  some  polymers  than  rock.  For  comparison, 
the  Poisson's  ratio  of  the  Esopus  shale  is  0.27.  Secondly, 
it  was  very  difficult  to  propagate  S  waves,  particularly  those 
in  the  direction  normal  to  the  bedding  (plane  of  polarization 
parallel  to  bedding.)  Furthermore,  a  frequency  dependent 
component  of  the  spectral  ratios  was  observed.  This  was  not 
true  for  the  Esopus  shale. 

Q  versus  confining  pressure  (Pf  =  0)  for  both  P  and  S 
waves  in  both  directions  are  shown  in  Figure  19.  The  results 
are  peculiar  -  not  all  expected.  First  note  that  Q  in  all 


cases  remain  unchanged  with  pressure  -  confirming  the  lack 
of  cracks.  Qs  is  very  low,  as  would  be  expected,  and  Qs 
for  the  perpendicular  direction  is  lower  than  the  parallel 
one,  mimicking  the  anisotropic  velocity  behavior.  However, 
the  anisotropy  of  Q  in  P  waves  is  reversed  from  that  of 
velocity.  Furthermore,  there  is  a  large  discrepancy  in  the 
differences  between  Qp  and  Qs  for  the  two  orientations.  In 
the  parallel  case  Qs  is  nearly  equal  to  Q  .  Finally,  this  is 
the  only  case  we  have  seen  in  which  the  changes  in  Q  with 
pressure  are  less  than  the  changes  in  velocity. 
r. .  Tuffs 

The  study  of  attenuation  in  tuffs  was  prompted  by 
interesting  differences  in  pore,  grain  size  and  crack 
distributions  compared  to  the  sandstones  implied  by  static 
measurements.  It  was  found  that  the  tuffs  contain  fewer 
small  aspect  ratio  cracks  than  either  the  Berea  or  Navajo 
sandstones  even  though  the  total  porosities  are  roughly 
equivalent.  Also,  the  observed  attenuation  in  both 
Ammonia  Tanks  and  Tunnel  Beds  tuffs  appear  to  be  dominated 
by  scattering  even  though  the  effect  of  cracks  is  evident. 


In  Tunnel  Beds,  the  Q  measurements  are  further  complicated 
by  the  presence  of  pore  collapse  under  pressure.  For  both 
tuffs,  only  dry  samples  were  examined. 

P  and  S  wave  velocities  for  the  Ammonia  Tanks  tuff  are 
shown  as  functions  of  confining  pressure  (Pf  =0)  in 
Figure  20.  Velocity  increases  more  gradually  with  pressure 
than  for  the  sandstones  and  does  not  begin  to  level  off  until 
about  1500  bars.  P  wave  velocities  are  similar  in  magnitude 
to  the  Navajo  and  about  5%  higher  than  the  Berea  even  though 
the  measured  porosity  is  much  less,  5%  compared  to  16-18%. 

S  wave  velocities  are  uniformly  lower  than  the  sandstones  by 
hbcut  5%.  This  implies  isolated  pores  and  cracks  mav 

exist  in  the  rock  and  that  the  actual  porosity  is  higher  than 
that  measured. 

The  Q  analysis  was  complicated  by  the  presence  of 
scattering  which  introduces  a  frequency  dependent  component 
to  the  apparent  attenuation.  Thus  the  spectral  ratios  are 
not  linear  over  a  broad  frequency  range.  The  Q  values 
calculated  are  strictly  applicable  only  to  a  frequency 
bandwidth  limited  to  the  peak  power  of  the  spectrum.  For 
the  Ammonia  Tanks,  this  occurs  at  about  0.4  MHz  for  P  waves 
and  0.5  MHz  for  S  waves.  The  resulting  Q  values  are  shown 
in  Figure  21.  The  most  outstanding  features  of  the  data 
are  the  low  Q  values  at  all  pressures  and  that  Qs  is  nearly 
1.5  times  Qp.  Since  one  would  intuitively  expect  that  the 
attenuation  due  to  scattering  would  be  relatively  constant 


with  pressure,  the  increase  in  Q  and  Q  with  pressure  may 

P  s 

be  due  to  crack  closure  just  as  with  the  sandstones.  However, 
as  with  velocity,  the  rate  of  increase  is  slow  compared  to  the 
sandstones. 

The  large  difference  between  Q  and  Q  (also  observed  in 

s  p 

the  Tunnel  Beds  as  will  be  seen)  is  obtainable  from  the  Walsh 
(1966)  friction  mechanism  but  requires  a  lower  coefficient 
of  friction  than  for  the  dry  sandstone  since  the  Poisson's 
ratio  is  0.24  compared  to  0.17.  Another  explanation  may  be 
that  the  attenuation  due  to  scattering  of  incident  S  waves  is 
less  than  for  incident  P  waves. 

T’he  Tunnel  Beds  tuff  was  more  difficult  to  work  with  compared 
to  the  Ammonia  Tanks.  With  its  very  high  porosity 
(30%) ,  vugular  pore  spaces,  and  large  grain  sizes,  scattering 
completely  dominates  the  anelastic  behavior  at  ultrasonic 
frequencies.  Pore  collapse,  initiated  at  about  800  bars  and 
observed  in  the  velocities  shown  in  Figure  22,  further 
complicates  the  picture.  P  and  S  wave  velocities,  about 
25-30%  lower  than  the  Ammonia  Tanks  do,  however,  behave 
predictably  and  in  a  manner  similar  to  the  Bedford  limestone. 

Under  loading,  velocity  increases  slightly  at  low  pressures 
due  to  crack  closure  and  then  levels  off.  At  800  bars,  velocity 
increases  again  as  pores  collapse.  During  unloading  of  the 
sample,  velocities,  are  higher  than  the  loading  cycle  for 
pressures  greater  than  800  bars  but  are  lower  at  lower 


pressures  due  to  cracks  generated  during  pore  collapse. 


Q  values  were  again  determined  for  a  limited  bandwidth 
centered  on  roughly  0.4  MHz  for  both  P  and  S  waves  and  are 
plotted  in  Figure  23.  Q  is  essentially  constant  with 
pressure,  consistent  with  the  domination  of  scattering. 

There  is,  however,  evidence  of  the  effect  of  pore  collapse 
barely  observable  in  the  loading  cycle  but  more  evident  in 
the  unloading  cycle.  A  clear  decrease  in  Q  relative  to  the 
loading  path  that  can  be  attributed  to  crack  production  during 
pore  collapse.  As  with  the  Bedford  limestone,  the  difference 
seems  to  be  greater  for  P  waves.  Also,  as  seen  for  the  Ammonia 


Tanks,  Q  is  less  than  Q_  by  about  the  same  amount.  Here, 
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CONCLUSIONS 

Ultrasonic  P  and  S  wave  Q  values  have  been  presented  for 
a  suite  of  dry  and  saturated  rocks  including  sandstones, 
a  limestone,  shales  and  tuffs  as  functions  of  pressure.  These 
measurements  showed  the  effects  of  different  rock  properties 
on  the  attenuation  of  seismic  waves .  Several  general  aspects 
concerning  the  behavior  of  attenuation  under  the  experimental 
conditions  studied  in  this  paper  are: 

1.  Q  increases  (attenuation  decreases)  as  a  function  of 
confining  pressure  (or  differential  pressure  for  saturated  rocks) . 
The  rate  of  increase  varies  depending  on  rock  type  and  crack 
density  and  distribution. 

2.  For  the  sandstones  and  limestone  measured  under  dry 

and  water  saturated  conditions,  both  Q  and  Q  are  higher 

P  s 


for  the  dry  case  as  compared  to  the  saturated  values. 

3.  For  the  dry  sandstones,  Qs  is  about  equal  to  or 
slightly  higher  than  Qp.  Qg  is  larger  than  Qp  by  nearly 
a  factor  of  two  in  the  dry  tuffs. 

4.  For  all  the  saturated  rocks  (including  the  shales), 

Qs  is  lower  than  Qp. 

5.  The  pressure  effects  are  dramatically  enhanced  during 
pore  collapse  as  evident  from  both  the  velocity  and  attenuation 
data  for  the  Bedford  limestone  and  Tunnel  Beds  tuff. 

6.  The  attenuation  in  tuffs  appears  to  be  controlled  by 

factors  other  than  those  affecting  the  attenuation  in  other 
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Q  in  the  tuffs,  associated  with  scattering  of  the  ultrasonic 
wave  points  out  the  fact  that  two  different  rock  types  with 
similar  porosities  may  exhibit  large  differences  in  attenuation 

7.  The  behavior  of  attenuation  in  shales  also  seems  to 
indicate  a  different  mechanism  than  might  be  proposed  for  the 
sandstones  or  limestone.  In  particular,  the  lack  of  change 
in  Q  with  pressure  for  the  oil  shale  suggests  that  cracks  are 
not  involved. 

8.  Anisotropy  in  velocity  is  probably  accompanied  by 
anisotropy  in  attenuation.  As  illustrated  by  the  oil  shale 
data  though,  the  anisotropic  behavior  of  attenuation  may 
not  correspond  directly  to  that  of  velocity. 


The  comparison  of  the  data  presented  in  this  paper  to 
in  situ  or  other  laboratory  measurements  of  attenuation 
must  be  approached  with  caution.  Ultrasonic  attenuation  data 
for  the  Berea  sandstone  (Tokstiz  et  al. ,  1979)  similar  to  the 
results  shown  here,  have  been  modelled  in  terms  of  several 

t 

attenuation  mechanisms  (Johnston  et  al. ,  1979)  including 
sliding  friction  at  cracks  and  grain  boundaries,  (Walsh, 

1966)  ,  Biot  (1956a, b;  1962a, b)  type  fluid  flow,  squirting 
flow  from  thin  cracks  to  larger  pores  (Mavko  and  Nur,  1975; 
O'Connell  and  Budiansky,  1977),  viscous  shear  relaxation  in 
fluid  saturated  cracks  and  pores  (Walsh,  1968  and  1969)  and 
scattering  from  pores  in  the  rock  (Yamakawa,  1962).  These 
theoretical  models  may  be  used  to  extrapolate  the  ultrasonic 
data  to  other  frequencies. 

While  the  frequency  dependence  of  attenuation  may 
possibly  be  determined,  one  aspect  of  the  problem  not 
considered  in  this  paper  is  the  variation  of  Q  with  strain 
amplitude.  Numerous  experimental  studies  have  shown  that 
although  Q  is  relatively  independent  of  amplitude  for  strains 
less  than  about  10~^  it  decreases  rapidly  at  higher  strains 
(Gordon  and  Davis,  1968;  Gordon  and  Radu,  1971;  Peselnick, 
and  Outerbridge,  1961;  Winkler  and  Nur,  1978).  Johnston 
(1978)  has  shown  that  apparent  discrepancies  in  Q  measurements 
obtained  using  low  amplitude  resonant  bar  techniques  and 
ultrasonic  methods  may  in  part  be  explained  by  amplitude 


dependence.  Thus,  if  amplitudes  at  seismic  frequencies  are 
different  from  those  used  in  laboratory  measurements,  absolute 
estimates  of  Q  are  invalidated. 

Further  work  is  required  to  characterize  the  differences 
in  Q  within  specific  rock  types ,  to  investigate  the  nature  of 
crack  lubrication  due  to  fluid  saturation  and  to  study  the 
effects  of  anisotropy  such  as  observed  in  the  oil  shale. 

These  could  eventually  lead  to  predictive  models  of 
attenuation  in  rocks  providing  information  on  both 
microstructure  and  fluid  content. 
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Table  1.  Sample  Parameters 


Berea  Sandstone 


Navajo  Sandstone  16.4 

Bedford  Limestone  11.9 
Esopus  Shale 
Colorado  Oil  Shale 
Ammonia  Tanks  Tuff  5 . 8 

Tunnel  Beds  Tuff  29.6 


Total  Porosity  Permeability  Grain  Pbulk  Pgrain 
%  md  Diameter  _  / _ 3 


18.4  75  0.10  2.137  2.620 

16.4  100  0.15  2.189  2.620 

11.9  1  0.75  2.335  2.652 

-  2.606  - 
»  -  -  2.004  - 

:  5.8  -  1-2  2.314  2.457 

29.6  -  2-5  1.607  2.283 
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FIGURE  CAPTIONS 


Fig.  la.  Attenuation  characteristics  of  P  waves  in  the  dry 
Nava- o  sandstone  at  two  confining  pressures,  Pc  =45 
and  300  bars  (Pf  =  0).  Left:  Pulse  shapes  in  the 
sandstone.  Middle :  Normalized  Fourier  amplitude  as  a 
function  of  frequency  for  the  rock  (solid  line)  and 
aluminum  standard  (dashed  line).  Right:  Natural 
logarithm  of  the  aluminum  (Al)  to  rock  (A2)  amplitude 
ratios  as  a  function  of  frequency.  Points  are  actual 
ratios  and  the  dashed  line  is  the  linear  fit  to  the  data. 
The  slope  of  this  line  is  proportional  to  1/Q. 

Fig.  lb.  Attenuation  charat-teristi  rs  of  S  waves  in  the  drv 
Navajo  sandstone.  See  Fig.  la  for  explanations. 

Fig.  2.  Linear  strain  as  a  function  of  hydrostatic  confining 
pressure  in  the  Berea  sandstone.  Loading  and  unloading 
paths  are  denoted  by  the  arrows . 

Fig.  3.  Linear  strain  as  a  function  of  hydrostatic  confining 
pressure  in  the  Bedford  limestone  for  two  successive 
pressure  cycles. 

Fig.  4.  Congressional  (P)  and  shear  (S)  wave  velocities  for 

the  dry  (open  symbols)  and  water  saturated  (solid  symbols) 
Berea  sandstone  as  a  function  of  differential  pressure. 

As  with  the  following  plots  of  velocity,  the  errors  are 
within  the  size  of  the  symbols  and  thus  no  error  bars  are 
drawn.  These  velocities  are  used  for  the  Q  calculations. 


Fig.  5.  Q  values  of  P  and  S  waves  as  a  function  of  confining 
pressure  in  the  dry  Berea  sandstone  (Pf  =0).  In  this 
and  later  plots  of  Q,  representative  error  bars  are  shown. 

Fig.  6.  Q  values  of  P  and  S  waves  as  a  function  of  differential 
pressure  in  the  saturated  Berea  sandstone  (Pf  =  0.465  Pc). 

Fig.  7.  Q  values  of  P  and  S  waves  of  a  function  of  confining 
pressure  in  the  partially  saturated  Berea  sandstone 
(Pf  =  0) . 

Fig.  8.  P  and  S  wave  velocities  as  a  function  of  fluid 
pressure  and  calculated  gas  saturation  for  the  gas 
release  experiment  in  the  Berea  sandstone.  The  dashed 
repres^n  theoretical  ly  calculated  P  wave 

velocity  for  a  fully  saturated  rock.  The  differential 
pressure  is  200  bars. 

Fig.  9.  Q  values  of  P  and  S  waves  as  a  function  of  fluid 

pressure  and  calculated  gas  saturation  in  the  Berea 

sandstone  (P.  =  200  bars), 
a 

Fig.  10.  P  and  S  wave  velocities  as  a  function  of  differential 
pressure  in  the  dry  and  saturated  Navajo  sandstone. 

Fig.  11.  Q  values  of  P  and  S  waves  as  a  function  of  confining 
pressure  in  the  dry  Navajo  sandstone  (Pf  =  0). 

Fig.  12.  Q  values  of  P  and  S  waves  as  a  function  of  differential 
pressure  in  the  saturated  Navajo  sandstone  (P^  =  0.465  Pc) . 

Fig.  13.  P  and  S  wave  velocities  as  a  function  of  differential 
pressure  in  the  dry  and  saturated  Bedford  limestone.  The 
unloading  path  is  denoted  by  the  dashed  lines. 


Fig.  14.  Q  values  of  P  and  S  waves  as  a  function  of  confining 
pressure  in  the  dry  Bedford  limestone  (Pf  =  0) .  The 
unloading  path  for  the  P  wave  is  denoted  by  the  dashed 
line. 

Fig.  15.  Q  values  of  P  and  S  waves  as  a  function  of  differential 
pressure  in  the  saturated  Bedford  limestone  (P^  =  0.465  Pc) . 
Unloading  paths  are  denoted  by  dashed  lines. 

Fig.  16.  P  and  S  wave  velocities  as  a  function  of  confining 
pressure  in  the  saturated  Esopus  shale  (Pf  =  0) .  The 
sample  was  held  at  1000  bars  for  16  hours  resulting  in  the 
velocity  increase  due  to  pore  pressure  equilibration, 
rsnshnd  lines  denote  unloading  path?. 

Fig.  17.  Q  values  of  P  and  S  waves  as  a  function  of  confining 
pressure  in  the  saturated  Esopus  shale  (Pf  =  0) .  Dashed 
lines  denote  unloading  paths. 

Fig.  18.  P  and  S  wave  velocities  as  a  function  of  confining 
pressure  in  the  Colorado  Oil  Shale  (Pf  =  0)  for  the 
directions  perpendicular  and  parallel  to  the  bedding  planes. 

Fig.  19.  Q  values  of  P  and  S  waves  as  a  function  of  confining 
pressure  in  the  Colorado  Oil  Shale  (P^  =  0)  for  directions 
perpendicular  and  parallel  to  the  bedding  planes. 

Fig.  20.  P  and  S  wave  velocities  as  a  function  of  confining 
pressure  for  the  dry  Ammonia  Tanks  tuff  (Pf  =  0) . 
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21.  Q  values  of  P  and  S  waves  as  a  function  of  confining 
pressure  in  the  dry  Ammonia  Tanks  tuff  (Pf  =  0). 

22.  P  and  S  wave  velocities  as  a  function  of  confining 
pressure  in  the  dry  Tunnel  Beds  tuff  (Pf  =  0) .  Dashed 
lines  denote  unloading  paths. 

23.  Q  values  of  P  and  S  waves  is  a  function  of  confining 
pressure  in  the  dry  Tunnel  Beds  tuff  (Pf  =  0) .  Dashed 
lines  denote  unloading  paths. 
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